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discovery that neutrinos are massive particles, that has implications from high-energy 
Q \ physics to cosmology. A new method for the production of intense and pure neutrino 

^ | ' beams has been proposed recently: the "beta-beam". It exploits boosted radioactive 

ions decaying through beta-decay. This novel concept has been the starting point 
| for a new possible future facility. Its main goal is to address the crucial issue of the 

existence of CP violation in the lepton sector. Here we review the status and the recent 
, developments with beta-beams. We discuss the original, the medium and high-energy 

' scenarios as well as mono-chromatic neutrino beams produced through ion electron- 

. capture. The issue of the degeneracies is mentioned. An overview of low energy 

beta-beams is also presented. These beams can be used to perform experiments of 
■ interest for nuclear structure, for the study of fundamental interactions and for nuclear 

\Q , astrophysics. 
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1. Introduction 

The observations made by the Super-Kamiokande (Fukuda et al 1998), the K2K 
(Ahn et al 2003), the SNO (Ahmad et al 2001) and the KAMLAND (Eguchi et 
al 2003) experiments have brought a breakthrough in the field of neutrino physics. 
The longstanding puzzles of the solar neutrino deficit, observed by Davis' pioneering 
measurement (Davis 1964) and of the atmospheric anomaly have been clarified. We 
know now that the expected fluxes are reduced due to the neutrino oscillation 
phenomenon: the change in flavour that neutrinos undergo while traveling, first proposed 
by Pontecorvo (Pontecorvo 1957). An impressive progress has been achieved in our 
knowledge of neutrino properties, thanks to these experimental advances. Most of the 
parameters of the Maki-Nakagawa-Sakata-Pontecorvo (MNSP) unitary matrix (Maki 
et al 1962), relating the neutrino flavor to the mass basis, are nowadays determined, 
except the third neutrino mixing angle, usually called #i 3 . However, this matrix might 
be complex, meaning there might be additional phases - one Dirac and two Majorana 
phases, if neutrinos are Majorana particles. A non-zero Dirac phase introduces a 
difference between neutrino and anti-neutrino oscillations and implies the breaking of 
the CP symmetry in the lepton sector. On the other hand, the CP violating Majorana 
phases can manifest themselves in other processes like e.g. neutrinoless double-beta 
decay. The Majorana or Dirac nature of neutrinos, the value of the neutrino magnetic 
moment, the existence of sterile neutrinos, the neutrino mass scale (and hierarchy) and 
the possible existence of CV violation are among the crucial issues that remain open. 

Neutrino oscillations imply that neutrinos are massive particles, contrary to what 
has been believed so far, and represent the first direct experimental evidence for physics 
beyond the Standard Model of elementary particles and interactions. Understanding 
the mechanism for generating the neutrino masses and their small values is clearly a 
fundamental question, that needs to be understood. On the other hand, the presently 
known (as well as unknown) neutrino properties have important implications for other 
domains of physics as well, among which astrophysics, e.g. for our comprehension of 
processes like the nucleosynthesis of heavy elements, and cosmology, for example for 
unraveling the asymmetry between matter and anti-matter in the Universe. Stars and 
the Universe constitute a precious laboratory for the study of the yet unknown neutrino 
properties, of fundamental interactions and for the search of new physics (Raffelt 1996). 
The present limits on the neutrino magnetic moment, on sterile neutrinos, or on the 
neutrino energy density coming from the cosmic microwave background and large scale 
structures (Bilenky et al 2003) are talkative examples of the strong constraints that can 
be extracted from the field of neutrino astronomy. 

If the study of astrophysical and cosmological phenomena offers a unique way to 
progress on the open questions, performing terrestrial - accelerator and non-accelerator - 
experiments represent an irreplaceable procedure to obtain direct (model-independent) 
information. Many experiments, running, planned or under study, will bring invaluable 
results in the near future. In particular, the measurement of the third, less known, 
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neutrino mixing angle, as well as of the Dirac phase need the availability of very intense 
neutrino beams, and might require new methods besides conventional beams (coming 
from the decay of pions and muons) pushed to their ultimate intensities, or super-beams. 
The idea of a neutrino factory was first put forward, which is based on the production, 
collection, acceleration and storage of muons to obtain muon and electron (anti)neutrino 
beams. The neutrino energies are of several tens of GeV and the needed baseline of the 
order of 3000 Km. The aim is to address in one step the leptonic CP study, the neutrino 
hierarchy and the precise measurement of the #23, Am^ and #13 oscillation parameters 
(Albright et al 2004). 

Recently Zucchelli proposed the idea of producing electron (anti)neutrino beams 
using the beta-decay of boosted radioactive ions: the "beta-beam" (Zucchelli 2002). 
Such a novel concept has three main advantages: well-known fluxes, purity (in flavour) 
and collimation. This simple and attractive idea exploits major developments in the 
field of nuclear physics, where radioactive ion beam facilities, either planned or under 
study, are expected to reach ion intensities of 10 11-13 per second. The main goals are to 
perform studies in the fields of nuclear structure, of fundamental interactions, in nuclear 
astrophysics, solid state physics, as well as for medical applications. Such intensities can 
render neutrino accelerator experiments using ions feasible as well, thanks to the beta- 
beam concept. 

In the original scenario proposed by Zucchelli, the ions are produced, collected 
and accelerated first at several hundred MeV per nucleon, then at several tens GeV 
energies - after injection in the Proton Synchrotron (PS) and Super Proton Synchrotron 
(SPS) accelerators^] - and finally stored in a storage ring of 7.5 km total length and 
of 2.5 km straight sections. This facility is based on reasonable extrapolation of 
existing technologies and exploits already existing accelerator infrastructure to reduce 
cost. The neutrino beam produced by the ions decaying along the straight sections 
point to a gigantic Cerenkov detector, located at the (upgraded) Frejus Underground 
Laboratory, in order to study CV violation, through a comparison of u e — *> and 
- ► V\i oscillations. A first feasibility study is performed in (Autin et al 2003) and 
discussed in (Lindroos 2003, Benedikt et al 2004). The discovery potential is now 
actively investigated (Mezzetto 2003, Bouchez et al 2003, Mezzetto 2005, Mezzetto 
2006, Burguet-Castell et al 2004, Burguet-Castell et al 2005, Campagne et al 2006) and 
is also discussed in (Guglielmi et al 2005). Other technologies are also being considered 
for the detector, like liquid Argon Time Projection Chambers (Rubbia C 1977, Rubbia 
A 2004, Ereditato and Rubbia 2004 and 2006) and iron calorimeters. 

The interest of the beta-beam concept for establishing a facility producing low- 
energy neutrino beams is soon recognized (Volpe 2004). Here the ions are accelerated at 
7 = 5 — 14 (7 is the Lorentz ion boost) to produce neutrino beams in the 100 MeV energy 
range. This idea opens new axis of research compared to the original scenario. Such 
beams can in fact be used for the study of nuclear structure, of fundamental interactions, 

X Zucchelli considers CERN as a possible site for a beta-beam facility (Zucchelli 2002). 
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and of nuclear astrophysics. In particular, the physics potential covers experiments on 
spin-isospin and isospin nuclear excitations and neutrino-nucleus interactions (Volpe 
2004, Serreau and Volpe 2004, McLaughlin 2004, Volpe 2005a), measurements of the 
Weinberg angle (Balantekin et al 2006a), of the neutrino magnetic moment (McLaughlin 
and Volpe 2004) and CVC tests (Balantekin et al 2006b), or core-collapse Supernova 
physics (Volpe 2004, Jachowitz and McLaughlin 2005 and 2006). A small devoted 
storage ring appears as more appropriate for these low energy applications (Serreau and 
Volpe 2004). Its feasibility study is now ongoing (Chance and Payet 2005). 

Other scenarios for the study of CV violation have been proposed afterwards, where 
the energy of the ions is much higher, the 7 ranging from 150 (Burguet-Castell et al 
2005) to several hundreds (medium energy) to thousands (high energy) (Burguet-Castell 
et al 2004). The value of 150 GeV per nucleon comes from the maximum acceleration 
that can be attained in the SPS. The baseline scenario in this case is the same as 
the original one. On the contrary, the medium and high energy options require major 
changes in the accelerator infrastructure, such as a refurbished SPS (or even the LHC) 
at CERN, as well as bigger storage rings. An increased distance between the source 
and the detector is also necessary to match the same oscillation frequencies. Therefore 
such scenarios imply further locations for the far detector, such as the Canfranc or 
the Gran Sasso Underground Laboratories. The physics potential is being actively 
investigated (Terranova et al 2004, Huber et al 2006, Agarwalla et al 2005, Donini et al 
2005a and 2006) and is discussed in (Terranova 2005, Gomez-Cadenas 2005, Migliozzi 
2005). It covers the third neutrino mixing angle, the CV violating phase, as well as the 
neutrino mass hierarchy. Some reduction of the degeneracy problem is also expected. 
This enlarged discovery potential is due to the information that can be extracted on 
the energy of the events; while in the original scenario the experiment is a counting 
measurement only. However, even if from the physics point of view pushing to higher 
energies appears as attractive, a specific and extensive feasibility study is needed, to 
determine in particular the ion intensities, that drastically influence the sensitivities, 
as well as the characteristics of the storage ring. From this point of view, the medium 
energy option appears as more promising, whereas the high energy one keeps at a much 
more speculative level. 

The presence of degeneracies constitute an important open issue in the analysis 
of future accelerator neutrino experiments. This consists in the appearance of several 
disconnected regions in the oscillation parameter space, besides the true solutions. The 
origin of such regions is both the inherent three flavour structure of the oscillation 
probabilities, and some still unknown neutrino properties, namely the sign of Am^ 
and the octant of #23- The "fake" regions are usually referred to in the literature as 
the sign, the octant and the mixed degeneracies (or clones). They might render the 
identification of the true values of 5 and #13 more complex, depending on our knowledge 
of these unknowns at the time a beta-beam experiment takes place. Several possible 
strategies can be followed to resolve the problem of the clone solutions. An extensive 
analysis in the context of beta-beams is performed in (Donini et al 2004, Donini et al 
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2005b, Donini et al 2005c, Campagne et al 2006, Donini and Fernandez-Martinez 2006, 
Fernandez-Martinez 2006). An interesting and promising possibility is pointed out in 
(Huber et al 2005b), i.e. to combine accelerator experiments with atmospheric data in 
the same detector. 

Finally, monochromatic neutrino beams produced by boosted ions decaying through 
electron capture are proposed in (Bernabeu et al 2005a) and further discussed in 
(Bernabeu et al 2005b, 2005c and 2006). The baseline envisaged here is the same 
as for the original beta-beam. In order to search for CV violation, instead of comparing 
neutrino versus anti-neutrino oscillations, a comparison of u e — > oscillations (only 
v e are available) at different neutrino energies is necessary. Again, a specific feasibility 
study is needed in particular to state the ion rates at production and in the storage ring, 
since such a configuration requires the acceleration and storage of not fully stripped ions. 

An attractive aspect of a beta-beam facility is that it offers strong synergies between 
different communities, to realize a multidisciplinary program covering hot open issues in 
various domains. In particular, the synergy between the nuclear and the neutrino physics 
communities has been quickly recognized. This has merged in a common feasibility 
study within the EURISOL Design Study, financed by the European Community in the 
FP6 Program (2005-2009), that is now ongoing. On the other hand the far detector 
can be used for other purposes, in particular, to measure neutrinos from a core-collapse 
Supernova explosion (in or outside our Galaxy), and to improve the present sensitivity 
on the proton decay. Finally, it is interesting to note that beta-beams can produce 
neutron beams as well, thanks to the beta-decay of boosted radioactive ions which emit 
delayed neutrons. This interesting possibility deserves further investigation. 

Clearly the beta-beam concept has a broad physics potential, which explains why 
this field has been developing very fast since its very beginning, bringing many new ideas 
and scenarios. In this review we summarize the work done and try to give an overall view 
of the field. The review is meant not only for experts, but also as an introduction for non- 
experts. For this reason, many technical details are skipped, unless strictly necessary. 
From now on, unless stated differently, with the terminology "beta-beams" we refer to 
the original scenario. In Section I we recall the phenomenon of neutrino oscillations, the 
recent discoveries in neutrino physics, as well as how neutrino oscillations can be used to 
address the breaking of the CV, T and CVT symmetries in the lepton sector. Section II 
presents the original beta-beam concept and baseline scenario, as proposed by Zucchelli, 
as well as the following sensitivity studies. In Section III low energy beta-beams and the 
physics potential discussed so far are presented. Section IV is devoted to the medium 
and high energy scenarios. The issue of the degeneracies is mentioned in Section V. 
Mono-chromatic neutrino beams based on electron capture are discussed in Section VI. 
Conclusions and perspectives are made in Section VII. 
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1.1. The neutrino oscillation phenomenon 

Let us briefly remind the basic ingredients of neutrino oscillations. This phenomenon 
arises if neutrinos have masses and there is mixing. The evolution equation of the mass 
eigenstates in vacuum is simply i(d/dt)\v m ) = H m \v m ). This becomes i(d/dt)\u a ) = 
UH m U'\i/ a ) in the flavor basis, where U is the unitary matrix which relates the weak 
v a (a = e, /i,r) to the mass z/j (i = 1,2,3) eigenstates. Each flavour eigenstate can be 
written as \u a ) = U*^) . 

In general, a unitary n x n matrix depends on n(n — l)/2 angles and n(n + l)/2 
phases. In the Dirac case, 2n — 1 phases can be removed by a proper re-phasing of the 
left-handed fields, leaving (n — l)(n — 2)/2 physical phases. Therefore, CV violation is 
only possible in the case of n > 3 generations. In the Majorana case there is less freedom 
to re-phase the fields, since the phases of the neutrino fields cannot be absorbed in the 
Majorana mass terms. Only n phases can be removed, leaving n(n — l)/2 physical 
phases. Out of these (n — l)(n — 2)/2 are the usual Dirac phases, while n — 1 are specific 
to the Majorana case and are called Majorana phases. The latter do not lead to any 
observable effects in neutrino oscillations and are not discussed further in this review. 

In the simple case of two families, this matrix reduces to the rotation matrix: 

v e \ _ ( cos 9 sin9 \ / v x \ , , 

i/„ ) ~ \ -sin^ costf ) \ u 2 ) U 

where the angle 9 is called the neutrino mixing angle, since it mixes the mass eigenstates 
Vx with mass m\ and v 2 with mass m 2 . 

If one assumes that a given source is producing a neutrino flux of given flavor \v a ) 
at t = x = 0, the neutrino state at a later time t is then \v a (t)) = U^e~ lEit \vj) , and 
evolves as an admixture of the mass eigenstates. The probability amplitude of finding 
the neutrino at the time t in a flavour state v$ is 

A(u a ^u p ;t) = U Pj e- iE ^ j . (2) 

The neutrino oscillation probability, i.e. the transformation probability of a flavour 
eigenstate neutrino \u a ) into another one \vp) is then 

P(v a -> vfit) = \A(u a -> u p -t)\ 2 = \U Pj e- iE ^U* aj \ 2 (3) 

Therefore the probability for detecting a state of a different flavor v@, at a given distance 
L from the source, is given, for relativistic particles, by: 

P{y a -> vp) w sin 2 29 sin 2 (1.27Am 2 ^) (4) 

E 

where L is in m (km) and E in MeV (GeV). The mixing angle 9 Eqs.flT]) and (j4j) 
determines the oscillation amplitude, whereas the difference of the square of the masses 
Am 2 = m\ — m 2 , gives the oscillation frequency (Figure [1]). The neutrino sources 
available (the sun, the atmosphere, reactors and accelerators) have different typical 
neutrino average energies, and distances between the source and the detector, and allow 
to probe various oscillation frequencies. Oscillation searches are called appearance 
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Figure 1. Schematic drawing of the v e to appearance oscillation probability in the 
case of two families Eq.@. The oscillation amplitude and frequency depend on two 
neutrino parameters: 9 and Am 2 = m 2 , — m\, where mi and mi are the masses of the 
mass eigenstates j^iand v% respectively. 



experiments, if the appearance of a new flavour is looked for, or disappearance 
experiments, if a reduction in the neutrino flux from the source is searched for. 

Neutrino propagation through matter can modify oscillations significantly. This 
is known as the Mikheyev-Smirnov-Wolfenstein (MSW) effect (Mikheev and Smirnov 
1986, Wolfenstein 1978). Indeed, neutrinos of all flavors couple to electrons, protons 
and neutrons in matter through neutral currents. In the case of v e scattering through 
charged- currents is also effective. This scattering induces a matter-induced potential 
acting on the electron neutrinos, which depends on the electron density n e : 

V e cc = V2G F n e . (5) 

The neutral current contributions from electrons and protons cancel, if the medium 
is electrically neutral, whereas the scattering on neutrons induces equal potentials on 
all neutrinos. Staying in the two-flavor formalism, the time evolution of the flavor 
eigenstates is as followg§|: 




^g-cos29 + V2G F n e ^fsin2# 
"sin 20 ^# cos 29 




(6) 

4B — - AE ■ COy ZO I V ^ J 

Note that in the three-flavor formalism, the matter-induced potential modifies both v e 
versus and to v T oscillations, since these are coupled. Matter effects can enhance 
neutrino oscillations when the MSW resonance condition is met: 

/\ 771, ^ 

——cos26 = V2G F n e . (7) 
2hj 

This means, that independently on the value of the vacuum mixing angle 6, mixing in 
matter is maximal if the MSW condition Eq.([7]) is satisfied. 

§ Terms which contribute equally to the diagonal give equal phases and do not influence neutrino 
oscillations. 
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In the case of three families, the neutrino flavour and mass eigenstates are related 
through 



(8) 



This matrix is called the Maki-Nakagawa-Sakata-Pontecorvo (MNSP) matrix (Maki et al 
1962), the analogue of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and is usually 
parametrized as: 
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(9) 



where c,- 



cos 8ij and 



sin 9ij, 9ij are the rotation angles (Figure [2]) and 6 is a CP 
Dirac violating phase. If neutrinos are Majorana particles, the matrix U is multiplied 
by an additional matrix having only non-zero diagonal terms, i.e. diag(l,e tai ,e tot2 ), 
including two extra CP violating phases. As far as the oscillation frequencies are 
concerned, since Amj = 0, only two frequencies are independent. Note that the 
oscillation phenomenon gives information on the Am- only, while the neutrino mass 
scale needs to be obtained from other measurements, in particular from the endpoint 
spectrum of nuclear beta-decay - like the tritium based experiment KATRIN (Angrik 
et al 2004) - or from the possible observation of the neutrinoless double-beta decay in 
nuclei, where the latter is also sensitive to the neutrino nature and to the Majorana CV 
violating phases (see e.g. Elliott and Engel 2004). 

In the last years much progress has been obtained in our knowledge of neutrino 
properties (for a detailed description see e.g. Bilenky et al 2003). Two Am 2 and 



mixing angles are now known, i.e. 



Am 2 12 



7.9 10~ 5 eV 2 , sin 2 fl 



12 



0.3, and 
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Figure 3. Sensitivity on sin 2 29\z as a function of time. For each experiment, its 
sensitivity is shown (solid line) , as well as the world sensitivity without the experiment 
(dashed line). The overall world sensitivity as a function of time is also shown 
(Guglielmi et al 2005). 

Am 2 ^ = 2.4 10~ 3 eV 2 , sin 2 # 23 = 0.5, often referred to as the solar and the atmospheric 
oscillation parameters respectively, extracted from the SNO (Ahmad et al 2001) and 
KAMLAND (Eguchi et al 2003) experiments on one hand, and from the Super- 
Kamiokande (Fukuda et al 1998) and K2K (Ahn et al 2003) experiments on the other 
hand. The sign of the solar oscillation parameter Am 2 2 is known, while the atmospheric 
one Am 2 ^ remains unknown. This introduces degeneracies in the determination of the 
other unknown parameters (Section V). Note that the LSND collaboration has found a 
positive oscillation signal at Am 2 = 1 eV 2 as well (Athanassopoulos et al 1998 and 1996). 
The running miniBOONE and the future BOONE experiments will confirm/refute these 
observations. Waiting for these results, in this review, we will not consider the third 
oscillation frequency and will stay within the three flavour formalism. 

The precise value of the #13 angle is not yet known, the present upper limit 
sin 2 26* 13 < 0.1 is extracted from a reactor experiment called CHOOZ (Apollonio et 
al 1999). Future experiments will improve the present sensitivities on sin 2 2^13 and 
search first in the 10 _1 — 10~ 2 range (Figure [3]). The 10~ 3 — 10~ 4 range will only be 
reached with beta-beams while values smaller than 10~ 4 can be explored with a neutrino 
factory. It is important to remind that the knowledge of #13 is crucial for CP violation 
searches, since if by any chance nature has decided that (9 13 is strictly zero, there cannot 
be any CP violation through the Dirac 5 phase. 

1.2. CP , T \CPT violation searches 

The study of the weak interaction and of neutrino properties has brought some of the 
building blocks of the Standard Model. Since 1964, we know that CP is one of the 
symmetries broken by the weak interaction, as first observed in K mesons (Christenson 
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et al 1964), and recently in B mesons. The CVT theorem ensures that under given 
conditions (Lorentz invariance, locality of the interactions) a modern quantum field 
theory description of particles and their interactions preserves the CVT symmetry. This 
is the case for renormalizable abelian theories. However some variants of string theories 
predict a breaking of Lorentz invariance and therefore of the CVT symmetry. Note that, 
if CVT is conserved, an obvious consequence of CV violation is the violation of the T 
symmetry. 

Neutrinos can be used to study these symmetries and their breaking through the 
comparison of neutrino oscillations. If CV is conserved, the oscillation probability 
between particles and the corresponding antiparticles coincide: 

CV : P{y a ^v$t)=P{p a ^>vp;t). (10) 

The action of the particle-antiparticle conjugation on the lepton mixing matrix U Eq.Q 
amounts to the change U — > U*, meaning that CV is only conserved if the mixing 
matrix U is real or can be made real by re-phasing of the lepton fields. The search for 
the possible existence of CV violation in the lepton sector and of a non-zero the Dirac 
phase 5 Eq.flUJ) consists in searching for possible differences between neutrino versus 
anti-neutrino oscillations. 

Time reversal interchanges the initial and the final states, so if T is conserved one 

has: 

T : P(u a ^u /3 ;t) = P(vp^u a ;t). (11) 

Finally, the oscillation probabilities are invariant under the combined action of CV and 
T, if the following equality holds: 

CVT : P{v a -> up- 1) = P(P -> V a - 1). (12) 

The study of these symmetries requires in particular appearance oscillation 
experiments, since disappearance experiments like P(u e — > u e ) are time-reversal 
invariant, and clearly insensitive to CV breaking effects. For example, by using Eqs.flBJ) 
and (JBHHD, the — > v e appearance probability can be rewritten as (Richter 2000) : 



At??, 2 L 

P(^ -> v e ) = 4c 2 3 s 2 3 s 2 3 sin 2 X 



2a 



l + -^-2-(l-2s 

13 



2 



Am? v 13J 



(13) 



2 AmlL Am?,L Am? 2 L 

+8c 13 Si2Si 3 S23(Ci2C23COS 6 - S12S13S23) COS — — Sin — — Sin — — 

^t, V ^tjy ^tjy 

2 . Am^L Am?,L Am? 2 L 
— 0C1 9Ci 2 C23Si2Si3S23sm o sin sin sin 

13 23 13 23 A ^ A ^ A ^ 

Attt,' 2 L 

+4s 2 2 C 2 3 (c 2 3 C23 + S 12 S 23 S 13 _ 2C12C23S12S23S13 COS 5) SU1 JJ 



4E V 



9 9 9 Ami! Am?,L aL , 9 , 
ScLstst cos ^— sin — (1 - 2sfo). 

12 13 23 AE 4E 4 E \ U) 



The first line contains terms driven by #13, the second and third lines contain CV 
even and odd terms respectively, while the forth line is driven by the solar parameters 
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Figure 4. Drawing of the P{v^ — > f e ) oscillations as a function of the distance, 
assuming a monochromatic beam of 1 GeV neutrinos, and matching the solar oscillation 
frequency for 5 — (left figure) and the atmospheric oscillation frequency for 5 = —n/2 
(right figure). In the latter case, the different terms contributing to the oscillation 
probability Eq. (fT3]) are shown (Guglielmi et al 2005). 



(Figure H]). Matter effects developed at the first order are included in the last line, where 
a[eV 2 ] = ±2V2G F n e E u = 7M0- 5 p[g/cm 3 ]E v [GeV\. The CV odd term as well as matter 
effects change sign by replacing neutrinos with anti-neutrinos. The occurrence of matter 
effects, at baselines of about 0(1000) km (Figured]), allows to extract information on 
the mass hierarchy, since their sign depends on the sign of Am|. These effects can 
be distinguished from the true CV violating ones, because they have a different energy 
dependence. 

Clearly, the knowledge of the Dirac 5 phase depends on the knowledge of all other 
oscillation parameters at the time an experiment searching for CV violation takes place. 
As mentioned above, a particularly important piece of information is the value of the 
angle #13. Note that if #13 7^ 0, it might be difficult from one single experiment to 
unambiguously extract information about #13 and S because of the correlations and 
the degeneracies (Section V). Anyway, pure, intense and collimated neutrino beams are 
needed to determine very small #13 values and a non-zero phase, such as those produced 
with beta-beams. 

2. The original beta-beam scenario 

2.1. The beta-beam concept and neutrino fluxes 

When a nucleus undergoes a beta-decay, it emits either an electron neutrino, or an 
electron anti-neutrino. In the rest frame of the nucleus, such neutrinos are emitted 
isotropically. However, if the ions are accelerated, the momentum parallel to the ion 
beam gets boosted, while the perpendicular one keeps unchanged. As a consequence, 
the neutrino beam gets collimated with the emittance of the neutrino flux inversely 
proportional to the Lorentz gamma factor. A neutrino beam produced with this 
procedure is what we call, a "beta-beam" (Zucchelli 2002). The neutrino flux is well 
known, since the ion intensity can be precisely determined and the flux is pure in flavour. 
Its expression can be determined analytically as follows (Serreau and Volpe 2004). 
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Figure 5. Comparison of the v e — » z/^ and the P e — * appearance oscillation 
probability, as a function of the baseline distance L between the source and the 
detector, evaluated at the first atmospheric oscillation maximum E/L — \Am,2 3 \/2n. 
From about L = 0(1000) km, matter effects become significant and can be exploited 
to extract the sign of Am^. The solid and dashed curves show the true CV violation 
effects for 6 = 0° and 90°. 



The decay rate of a nucleus in the rest (cm) frame can be written as: 



dW 



dt 



<$> cm {E v )dE l 



d 2 Q 
An 



(14) 



where E v and Q denote respectively the energy and the solid angle of the emitted 
(anti)neutrino, and where the neutrino flux <& cm (E v ) is given by the well-known formula: 

<$> cm {E v ) = b El E e y/E*-m? e F{±Z, E e ) Q(E e - m e ) . (15) 

where the constant b = \n2/mlfti/ 2 , with m e the electron mass and ftyi the ft-value. 
The quantities appearing in the above expression are the energy E e = Q@ — E v of the 
emitted lepton (electron or positron), Qp is the Qp-value of the decay. The Fermi 
function F(±Z, E e ) accounts for the Coulomb modification of the spectrum. 

In the laboratory frame, where the boosted nucleus has a velocity v = f3c, the decay 
rate reads: 

dW 

lab 



dt 



1 d 2 Q 
-<$> lah (E u ,6)dE u °— 

7 4:71 



(16) 



where 7 = l/y/l — j3 2 is the time dilation factor and where E u and Q = (9,(p) now 
denote the energy and solid angle of the emitted (anti) neutrino in the laboratory (lab) 
frame, 6 being the angle of emission with respect to the beam axis. The boosted flux 
$ lab (E„,6) is given by: 



<S>lab(E u , 9) 



<P cm (E,7[l-/?cosfl]) 
7[l-/3cos0] 



(17) 
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The ions are stored in a storage ring of total length S with straight sections of length 
D. In the stationary regime the mean number of ions in the storage ring is 'yrg, where 
t = t 1 / 2 /ln2 is the lifetime of the parent nuclei and g is the number of injected ions 
per unit time. The total number of neutrinos emitted per unit time from a portion di 
of the decay ring is therefore 

dN v dW d£ 

— = 7tox — x— . (18) 
dt dt S K J 

Let us consider, for example, a cylindrical detector of radius R and depth h, aligned 

with one of the the storage ring straight sections, and placed at a distance d from the 

latter. After integration over the useful decay path and over the volume of the detector, 

the total number of events per unit time is: 

dN r°° 

= grnh x / dE v $ m {E v ) a{E v ) , (19) 
dt Jo 

where n is the number of target nuclei per unit volume, a(E v ) is the relevant interaction 
cross section with the target, and the total neutrino flux is 

$tot(E v ) = J — I — I -^^$ lab (E v ,9), (20) 

with 



S Jo h 



4 = 3+7+7' (21) 




D D.2 0.4 0.6 0.S 1 0.2 OA 0.6 0.S 1 



E, (GeV) (GeV) 

Figure 6. Comparison of neutrino fluxes from a super-beam (SPL) and a beta- 
beam. The neutrino beams are produced at CERN and sent to the Frejus Underground 
Laboratory, 130 km from CERN. Two options for the beta-beam are shown here. Left: 
The ions circulate together in the storage ring, with 7 = 60 (100) for 6 He ( 18 Ne) 
(Mczzetto 2005). Right: The ions circulate at the same 7 = 100, independently, in 
the storage ring. Note that the average neutrino energies are related to the ion boost 
through E v w 2-/Q p (Guglielmi et al 2005). 
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Figures [6] and [7] show examples of the neutrino fluxes from beta-beams where the 
ions are boosted at 7 of several hundreds or from low energy beta-beams with 7 = 7—14, 
respectively. 

2.2. The baseline 

The original conceptual design for a beta-beam beam facility includes three main steps 
(Zucchelli 2002, Autin et a/2002, Bouchez et al 2003, Benedikt et al 2004, Lindroos and 
Volpe 2004, Benedikt 2005). In the first stage protons delivered by a Superconducting 
Proton Linac (SPL) impinge on a target to produce radioactive species of interest, using 
an ISOL target. Such an intense proton driver (at present under study) would deliver 
2mA of 2.2 GeV (kinetic energy) protons, while the ISOL target would need 100 /j,A, 
that is 5 % of the proton intensity only. After ionization and bunching, the ions are 
accelerated to low energy, i.e. 300 MeV/A. This part of the facility is very close to the 
conceptual design of the EURISOL project (see EURISOL), aiming to the production 
and acceleration of very intense radioactive ion beams through the ISOL technique. In 
the second stage energies of about 100 GeV/A are attained using e.g. CERN accelerator 
infrastructures, like the PS and SPS. Finally the ions are stacked longitudinally in a 
new decay ring. Each of these steps require specific feasibility studies, in particular the 
bunching at low energy, space-charge limitations in the PS and SPS (keeping CERN 
in mind as a possible site), decay losses in the accelerator chain and, at high energy, 
the stacking process in the decay ring. Both the decay ring, the proton source and the 
low energy accelerators do not exist at CERN today. Nevertheless, the use of existing 
accelerator infrastructure represents an important saving. The beta-beam complex is 
shown schematically in Figure [H 

The choice of the isotopes is based on several criteria. The half-life should be such 
as to minimize the decay losses in the accelerator chain. Another constraint comes from 
the required collimation and energy of the neutrinos at the detector, which depends on 
the energy of the decaying ions through E u « These (and other) considerations 

point at present at two isotopes of particular interest, namely 6 He and 18 Ne as anti- 
neutrino and neutrino emitters respectively. Their lifetimes are 0.8 s and 1.6 s, while 
their Qp- values are of about 4 MeV. Both these ions can be produced in large quantities 
by the so-called ISOL method. These two noble gases have the advantage that they can 
easily diffuse from the target. For helium, it consists either of a water cooled tungsten 
core or of a liquid lead core which works as a proton to neutron converter surrounded 
by a beryllium oxide, aiming for 10 15 fissions per second. This yields a quasi-continuous 
beam, so the ions have to be bunched prior to further acceleration. Efficient bunching 
and full stripping of a high-intensity beam can be achieved using a high frequency ECR 
source (Sortais 2003). While such system does not exist today, theoretical calculations 
show that it can be constructed. Note that very intense 8 Li and 8 B beams can be 
achieved with a new "ionisation cooling technique" (Rubbia et al 2006). 

16 bunches (consisting for example of 2.5 10 12 ions each for helium) are first 
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Figure 7. Anti-neutrino fluxes from low energy beta-beams. The curves correspond 
to 7 — 7 (dotted line), 7 = 10 (broken line) and 7 = 12 (dash-dotted line) for the 
boosted hclium-6 ions. Similar curves are obtained for neutrinos from neon-18 ions. 
The detector is located at 10 m from a small devoted storage ring (see Section III) 
(Balantckin et al 2006a). 



accumulated in the PS, reduced to 8 bunches during the acceleration to intermediate 
energies, injected in the SPS and then in the storage ring in batches of four 10 ns bunches. 
This time structure is required to efficiently reduce the atmospheric background. Note 
that the SPS can accelerate 6 He and 18 Ne ions at a maximum 7 value of 150 and 
of 250 respectively (Lindroos 2005). Due to the relativistic time dilatation, the ion 
lifetimes reach several minutes, so that stacking the ions in the decay ring is mandatory 
to get enough decays and hence high neutrino fluxes. The ring has the shape of an 
hippodrome, with a total length of 6880 m and straight sections of 2500 m each (36 %). 
The average number of ions in the storage ring between two stacking cycles is given by 
N average = NbunchJT/T where Nf, unc h is the number of particles with mean- life r having 
Lorentz boost 7, injected into the ring every T seconds. 

The ion injection and merging with existing high density bunches in the decay ring 
represent an important challenge and require a new scheme. Stacking means a regular 
top-up of the stored beam before the existing particles cease to be useful. Many stacking 
methods involve some form of a beam cooling to increase phase space density and make 
room for more particles within a given transverse and longitudinal acceptance. For light 
ions at high energy, the classical methods of electron and stochastic cooling are excluded. 
Instead, a new scheme has been proposed combining bunch rotation and asymmetric 
bunch pair merging. The new bunches are off momentum and are injected in a high 
dispersion region on a matched dispersion trajectory. Subsequently, each injected bunch 
rotates a quarter turn in longitudinal phase space until the initial conditions for bunch 
pair merging are met. In the final step, each small fresh bunch is moved into the core of 
a large stored one and then phase space mixing occurs. The fact that only the central 
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Figure 8. Schematic drawing of a beta-beam facility complex at CERN. The 
production, collection and acceleration part shown on the left is very close to the 
EURISOL project (see EURISOL), at present under study, which aims at producing 
very intense radioactive ion beams. The ions are then injected in the Proton- 
Synchrotron (PS) and Super Proton-Synchrotron (SPS) accelerator infrastructures 
(central part) already existing at CERN, to achieve about a hundred GeV per nucleon 
(7 = 100). Finally the ions are stored to decay in an hippodrome- like storage ring 
(right part), that needs to be built. One of the long straight sections points to a far 
detector (Zucchelli 2002, Autin et al 2003). 

part of the existing bunch is combined with an incoming dense one results in a net 
increase in the core intensity of the resultant stored beam. The surrounding older ions 
are pushed out toward the bucket separatrix, where eventually the oldest are lost. This 
new method, called asymmetric bunch pair merging, has recently been demonstrated at 
the CERN PS. 

A principal difference between the acceleration to high energies of stable and 
radioactive ions is the additional losses due to the radioactive decays. The predicted 
intensities for 6 He and 18 Ne are shown in Table HJ The conventional losses through the 
accelerator chain are expected to be less than 50 % (based on operational experience at 
CERN). The isotopes proposed for the beta-beam conceptual design have been chosen 
such that no long-lived activity is left to contaminate the low energy machines. A first 
simulation of losses in the decay ring yields a dose rate in the arcs of 2.5 mSv/h after 30 
days of operation and one day of cooldown (Magistris and Silari 2003). Furthermore, 
the induced radioactivity in the ground water will have no impact on public safety. It 
has also been confirmed that the total loss in the accelerator chain will be beyond the 
nominal 1 W/m that permits hands-on maintenance (except for the PS). 
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Table 1. Intensities for the 6 He and 18 Ne ion beams along the acceleration chain. 
Only beta-decay losses are taken into account. The number of ions in the decay ring, 
including conventional losses, is 0.5N tot (Autin et al 2002, Lindroos 2003, Benedikt et 
al 2004, Benedikt 2005). 



Machine 


6 He ions extracted 


18 Ne ions extracted 


Batches 


Source 


2xl0 13 


8xlO n ions/s 


dc 


Rapid Cycling Synchrotron 


l.OOxlO 12 


4.09xl0 10 


16 


PS 


l.OlxlO 13 


5.19xlO n 


1 


SPS 


0.95xl0 13 


4.90xlO n 


OG 


Decay Ring (N tot ) 


2.02 xlO 14 


9.11X10 12 





2. 3. Sensitivity to CV violation 

The major goal of the beta-beam facility proposed by (Zucchelli 2002) is to measure 
CV violation in the lepton sector. This requires a comparison of neutrino versus anti- 
neutrino oscillations in an appearance experiment Eq. ffTOl) (note that disappearance 
experiments like v e — > v e are not sensitive to the CV violating phase 5). Since beta- 
beams produce pure electron neutrino (or anti-neutrino) beams, such experiment should 
search for differences in u e — *> versus v e — > oscillations Eq. (1131) . The scenario 
considered is to produce the neutrino beams at CERN (Figure [H]) and fire them to a 
water Cerenkov detector of 440 kton fiducial volume - like MEMPHYS (Mosca 2005a, 
Campagne et al 2006), UNO (Jung 2000), or HYPER-K (Itow et al 2001) - located 
at the (upgraded) Frejus Underground Laboratory, 130 km from CERN. If Lorentz ion 
boosts around 60-100 are chosen (FigJBJ), this distance magically fits the v e — > first 
oscillation maximum at the atmospheric scale (Figf4J). The detector is of the same kind 
of Super-Kamiokande (Fukuda et al 1998), but has a fiducial mass 20 times larger, 
requiring a cavern of about 1 million cubic meters. The ion intensity considered is 
2.9xl0 18 6 He and 1.1 xlO 18 18 Ne per year (Autin et al 2003). 

The physics potential of the original baseline scenario is being actively investigated 
(Mezzetto 2003, Bouchez et al 2003, Mezzetto 2005, Mezzetto 2006, Campagne et al 
2006) and is discussed in (Guglielmi et al 2005, Mosca 2005b, Bouchez 2005). The 
signal corresponds to the muons produced by charged- current events in water, mainly 
via quasi-elastic interactions at these energies. Such events are selected by requiring a 
single-ring event, with the same identification algorithms used by the Super-Kamiokande 
experiment (Fukuda et al 1998), and by the detection of the electron from the muon 
decay. At such energies the energy resolution is very poor due to the Fermi motion and 
other nuclear effects. For these reasons, a CV violation search with 7 = 60 — 100 is based 
on a counting experiment only. Energy cuts can be exploited when going to higher 7, 
as first suggested in (Burguet-Castell 2004) (see Section IV) e.g. to resolve some of the 
degeneracies and get information on the neutrino hierarchy. 

The beta-beam has no intrinsic backgrounds, contrary to conventional sources. 
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However, inefficiencies in particle identification, such as single-pion production in 
neutral-current v e (u e ) interactions, electrons (positrons) misidentified as muons, as well 
as external sources, like atmospheric neutrino interactions, can produce backgrounds. 
The background coming from single pion production has a threshold at about 450 
MeV, therefore giving no contribution for 7 < 55. Standard algorithms for particle 
identification in water Cerenkov detectors are quite efficient in suppressing the fake 
signal coming from electrons (positrons) misidentified as muons, especially if the signal 
for the muon decaying into electrons is also required. Concerning the atmospheric 
neutrino interactions, estimated to be of about 50/kton/yr, this important background 
is reduced to 1 event /440 kton/yr by requiring a time bunch length for the ions of 10 
ns. 

The choice of the ion acceleration is a compromise between having the 7 factor 
as high as possible, to profit of larger cross sections and better focusing of the beam 
on one hand, and keeping it as low as possible to minimize the single pion background 
and better match the CP odd terms Eq. (fl3i) on the other hand. If one requires the 
ions to circulate together in the decay ring (Figure ED, this constraints the ratio of 
the ion boosts to be 7( 6 He)/7( 18 Ne)=3/5. In this case the optimal boosts turn out to 
be 7( 6 He)=60 and 7( 18 Ne)=100, giving 0.24 GeV and 0.36 GeV as v e and v e average 
energies respectively (Figure [6] left). This choice is assumed in (Bouchez et al 2003, 
Mezzetto 2005). 

Table [2] shows the expected events with and without oscillations. One can see 
that the beam related background is zero, while some background events are expected 
especially for the case of neon-18, due to the single-pion misidentified events. The 
expected sensitivity on #13 is shown in Figure [9] after 5 years run, in comparison with 
planned or proposed experiments. Note that the value of this angle can be explored both 
with v e and v e disappearance and with and appearance experiments. Clearly, a 
beta-beam alone as well as a combination of a super-beam and a beta-beam sent to the 
same detector are powerful and allow to reach sensitivities on sin 2 2# 13 in the range of 
10 -3 — 10 -4 , improving the present CHOOZ upper limit (Apollonio et al 1999) by a 
factor of about 200. 

The sensitivity in the (#13, S) plane is shown in Figure dOj having fixed the 
other oscillation parameters to the values sin 2 2^23 = 1, Am^ = 2.5 10~ 3 eV 2 , 
sign(Am2 3 )=+l, sin 2 2# 12 = 0.8, Am 2 2 = 7.1 10 _5 eV 2 . The assumed corresponding 
errors are of 5 % and of 1% on Am| 3 and on sin 2 2#23 respectively, as expected from 
the T2K experiment (Itow et al 2001), while 10 % error is taken for the solar oscillation 
parameters. The curves are determined by fitting the number of muon-like events in 
the and oscillation experiments simultaneously (Table [2]). Correlations between 
#13 and 5 are fully accounted for and have a negligible effect (Figure [I0J see section V 
for a discussion on the intrinsic degeneracy). The ambiguity on the sign of Am| 3 is 
translated in a direct ambiguity on the sign of 5, since matter effects at 130 km are 
negligible Eq.t flBl) (Figure [5]). The ambiguity coming from the choice #23 or 7r/2 — #23 is 
formally taken into account but has no effects in the results shown here, since the choice 
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Table 2. Number of events expected after 10 years, for a beta-beam produced 
at CERN (Figure [5} and sent to a 440 kton water Cerenkov detector located at 
an (upgraded) Frejus Underground Laboratory, at 130 km distance. The results 
correspond to v e (left) and v e (right). The different 7 values are chosen to make 
the ions circulate together in the ring (Mezzetto 2005). 





6 He 


18 Ne 




(7 = 60) 


(7 = 100) 


CC events (no oscillation) 
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Oscillated (sin 2 29 13 = 0.12, 5 = 0) 
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Figure 9. Sensitivity on #13 (90 % C.L.) as a function of 6 for a beta-beam, a super- 
beam (SPL), and their combination. The present CHOOZ limit (Apollonio et al 1999) 
is shown, as well as the sensitivities of CNGS and JPARC (Migliozzi and Terranova 
2003) and BNL (Diwan et al 2003). The appearance results are calculated considering 
5 years measuring time (Mezzetto 2005). 



623 = 45° is made. A detailed discussion of the issues of degeneracies, correlations and 
clones is made in Section V. 

In order to achieve the highest sensitivity it is clear that the systematic errors 
have to be kept at their lowest possible level. Such errors come essentially from the 
knowledge of the cross sections - which at the neutrino energies of the beta-beams 
(several hundred MeV) are not well known - as well as from the precision with which 
all the other oscillation parameters are known, at the time such an experiment takes 
place. By placing a close detector of 1 kton at least at 1 km from the decay tunnel, one 
can exploit the same neutrino beams to perform precision cross section measurements 
(Mezzetto 2003). Note that, at the moment, the cross section uncertainty introduces 
variations in the sensitivities found in the literature (see e.g. Singh et al 2006a for a 
recent calculation of the neutrino-oxygen cross sections with beta-beams). As far as 
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Figure 10. Fits to #13 and S after ten years measurement time, computed with the 
events shown in Table [2j The stars indicate the parameters' values, while the different 
curves show the various confidence levels (Mezzetto 2005). 



the energy and the flux of the neutrino beams are concerned, one should keep in mind 
that contrary to conventional beams, in the case of beta-beams these are completely 
defined by the intensity of the ions and the kinematics Eqs. (fT6l) and ( 1TT1) . The flux at 
the far detector is also completely determined by the divergence of the beam which is 
perfectly known. One can therefore hope to measure the signal and the backgrounds 
with a residual systematic error of 2 %. Figure fTTl shows the importance of keeping the 
systematic errors at their lowest level. 



90 




0123456789 10 

ei3(cleg) 

Figure 11. Discovery potential for 5 at 3 a, considering 10 years measurement time, 
with the two ions running at 7 = 60, 100. The effect of systematic errors is shown. The 
curves correspond to 2 % (solid), 5 % (dashed), 10 % (dot-dashed) systematic errors. 
The results for a beta-beam with the two ions running at the same 7 = 75 with 2 % 
systematic error is also shown (dotted). The two grey curves present the sensitivity 
of a super-beam (solid) and super-beam combined to a beta-beam (dashed) with 2 % 
systematic error (Mezzetto 2005). 



A different optimization of the ion acceleration can be obtained if each type of ion 
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(either the v e or the v e emitter) circulates separately at its optimal 7 value, as suggested 
in (Mezzetto 2003, Burguet-Castell et al 2005, Lindroos 2005). In this case the single 
fluxes can be doubled by filling all the batches with the same type of ions (instead of 
alternating), which gives the same integrated number of decays per year. The physics 
potential obtained following such a prescription is studied with 7( 6 He) = 7( 18 Ne)=75 
in (Mezzetto 2003) and with 7 = 100 (Figure ED in (Mezzetto 2006, Guglielmi et al 
2005, Campagne et al 2006). These two choices well match the CERN-Frejus distance 
of 130 km. The sensitivity on the third neutrino mixing angle as well as on the CV 
violating phase is presented in Figures [12] and [13] respectively, showing an increased 
discovery potential for both observables compared to the case where the two types of 
ions circulate together in the storage ring. The beta-beam option alone is clearly a better 
strategy than a super-beam alone, as already apparent in Figures l9l and [TT1 Finally, the 
comparisons show that the largest discovery potential can be attained by a combination 
of a super-beam and a beta-beam. 

et 
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Figure 12. Sensitivity on the third neutrino mixing angle (90 % C.L.) as a function 
of the CP violating phase, with 2 % systematic errors, as in Figure [HI The discovery 
potential of Double CHOOZ (Ardellier et al 2004) is also presented. Here two different 
beta-beam options are compared, i.e. /3B60,ioo and /3Bioo,ioo where 7( 6 He)=60 and 
7 ( 18 Ne)=100 and 7( 6 He) = 7( 18 Ne)=100 are used respectively (FigureE). The results 
correspond to 5 years running time for the super-beam case (SPL-SB) with either 
2.2 or 3.5 proton kinetic energy, and 5 years v e and v e measurement time for the 
beta-beam case (Guglielmi et al 2005). 

Since a beta-beam would exploit at most 10 % of the SPL protons (the injector), 
one could produce both a beta-beam and a super-beam and send these neutrino beams 
to the same detector. This choice has several advantages. The use of the two beams does 
not only increase the statistics, but also gives the possibility to reduce the systematic 
errors and offers the necessary redundancy to firmly establish CV violating effects, within 
the reach of the experiments. Besides, the combination of a beta-beam and of a super- 
beam offers the largest discovery potential both with respect to the third mixing angle 
and to the CV violating phase. The last important advantage is that also T and CVT 
violation studies can be performed through a comparison of v e — > oscillations (from 
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beta-beams) to z/ M — > v e Eq. (|TTT) and to — > z/ e Eq.( fT2l) oscillations (from super-beams) 
respectively (Mezzetto 2003). 

An extensive study of the beta-beam option at 7 = 100 is made in (Campagne 
et al 2006) based on the GLoBES software (Huber et al 2005a), including correlations 
and degeneracies (Section V) and using atmospheric data in the analysis (Huber et al 
2005b). The GloBES software is a convenient tool to simulate long-baseline experiments 
and compare different facilities on the same footing. Correlations and degeneracies are 
fully taken into account letting all oscillations parameters to vary in the fits. A detailed 
comparison is made among the beta-beam, the super-bean[jj] (with 3.5 GeV proton 
kinetic energy), and T2K phase II (T2HK) facilities. The latter corresponds to the 
second phase of the T2K experiment in Japan, with a 4 MW proton driver and the 
Hyper-K detector located at 295 km (Itow et al 2001). It is a competing proposal with 
beta-beams with a similar size and timescale. In this work it is pointed out that if 
one assumes the CVT symmetry, a combination of a beta-beam and a super-beam sent 
at the same detector can be used to replace the anti-neutrino measurement by the T 
conjugate one. This implies that if both beams are available, the full information can 
be extracted with the neutrino beams only, reducing the 10 years measurement time 
down to 5. The corresponding sensitivity comes out to be very close to the 10 years 
measurement time with neutrinos and anti-neutrinos from one single experiment. 

Figures [TH and [15] show the sensitivity on the third mixing angle and the discovery 
reach with respect to CP respectively. One can see that the three facilities have quite 
similar performances and can clearly reach values as small as sin 2 2#i3 = 5 10 -3 , 
independently on 5. Besides, for a large fraction of 5 values, a sensitivity of 10 -3 
can be attained. Assuming 2 % systematical error, in the best case 5(7) 10~ 3 can be 

|| Note that the super-beam here has a new optimization of the beam energy as well as of the secondary 
particle focusing and decay. This new optimization has a significant impact on the physics performance 
(Campagne et al 2006). 



CONTENTS 



24 



7-10' 



3a discovery of a non-zero 9 

~l — i — i — i — i — | — i — i — i — 



13 



<x> 



1 -10 



3-10 




71 

true 5 



CP 



0.25 0.5 
fraction of true 5 



CP 



0.75 
values 



Figure 14. Sensitivity on #13 at 3cr for the (7 = 100) beta-beam (/3B), a super-beam 
(SPL), and T2HK (see text) as a function of the true value of i5cp (left) and as a 
function of the fraction of all possible values of Sep (right). The width of the bands 
corresponds to values with 2% to 5% systematical errors. The dashed curves correspond 
to the combination of the beta-beam and the super-beam with 10 yrs measurement 
time each and 2% systematical error (Campagne et al 2006). 



reached for the beta-beam and T2HK (SPL). In a combined beta-beam and super-beam 
(SPL) experiment with 10 years measurement time 4 10 -4 is reached for 45 % of 5 
values. Concerning CV, the maximal violation corresponding to 5 = tc/2 and 3ir/2 can 
be discovered down to sin 2 2# 13 = 6(8) 10~ 4 (99 % C.L.) for the beta-beam (super-beam); 
while the best sensitivity is obtained for sin 2 2^13 > 10~ 2 since at this value CV violation 
can be established for 73 % (75 %) for all 5 values with a beta-beam (super-beam). One 
can also see that the impact of systematical errors is crucial in the case of the T2HK 
experiment while the beta-beam and the SPL are much less sensitive. The comparison 
with T2HK shows that the 7 = 100 beta-beam option has a very similar discovery reach 
and offers a competitive strategy for future neutrino precision oscillation measurements. 

Most of the literature concerning beta-beams considers the water Cerenkov 
technique for the far detector. There are nowadays three proposal for very large water 
Cerenkov detectors: MEMPHYS in Europe (Mosca 2005a, Campagne et al 2006), UNO 
in the United States (Jung 2000) and HYPER-K in Japan (Itow et al 2001). This mature 
technology has already demonstrated excellent performances at the neutrino energies of 
several hundred MeV to 1 GeV. Other techniques are now being discussed, in particular 
the liquid Argon Time Projection Chamber (Rubbia C 1977, Rubbia A 2004, Ereditato 
and Rubbia 2004 and 2006) and iron calorimeters. It is important to emphasize that 
the large scale detector considered here can be exploited for a multidisciplinary program 
covering crucial issues in astrophysics, such as the detection of neutrinos from a(n) 
(extra-)galactic core-collapse Supernova, and in high energy physics, like proton decay. 
The field of neutrino astronomy has started with Davis' pioneering experiment (Davis 




Figure 15. Same as Figure [UJ but for the CV discovery reach. For parameter values 
inside the ellipse-shaped curves CP conserving values of 5cp can be excluded at 3a 
(Ax 2 > 9) (Campagne et al 2006). 



1964) which has measured solar neutrinos from the first time. Afterwards, neutrinos 
from a core-collapse Supernova have been observed during the explosion of the SN1987A, 
located in the Large Magellanic Cloud, at about 50 kpc from the Earth. For these 
important observations, both R. Davis and M. Koshiba received (with R. Giacconi) the 
Nobel Prize in 2002. Nowadays, neutrinos are currently used as a probe of astrophysical 
objects. In particular, the measurement of the neutrino lightcurve from a future core- 
collapse Supernova explosion constitute an essential ingredient to unravel the explosion 
mechanism, which keeps being a mystery: in the most recent tri-dimensional models 
the shock wave stalls and fails to eject the mantle. The neutrinos closely follow the 
explosion, from the collapse and the accretion phases, to the cooling of a proton-neutron 
star, or the formation of a black hole. Besides a recently proposed technique, that 
consists in adding Gadolinium to a water Cerenkov detector (Beacom and Vagins 2004), 
can increase the sensitivity to the diffuse Supernova neutrino background from past 
Supernova explosions, for which there is at present an upper limit. The observation of 
both future and past Supernova neutrinos would bring invaluable information not only 
for our knowledge of Supernova physics, but also on neutrino properties which are still 
poorly known (Raffelt 1996, Fogli et al 2005). The water Cerenkov detector considered 
for the beta-beam project can measure about 150000 events from a Supernova at 10 
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kpc - while presently running detectors can observe at most 10000 events - and about 
50 events from Andromeda. A flux of 250 events/10 years/500 kton from the diffuse 
background that can be increased by a factor of 10 with Gadolinium. Concerning proton 
decay one can attain a sensitivity up to 10 35 years in the e + 7r channel and a few 10 34 
years in the K + v channel (Jung 2000, Itow et al 2001, Mosca 2005a). 

3. Low energy beta-beams 

Very soon the idea of exploiting the beta-beam method to establish a facility that 
produces low energy neutrino beams - in the 100 MeV energy range - has been 
proposed, opening new axis of research in various domains of physics (Volpe 2004). 
This corresponds to a range of about 7-14 for the ion boosts 7 (Figure [7]). The physics 
potential of such a facility is now being analyzed. Several applications have been 
studied concerning nuclear structure studies and neutrino-nucleus interactions (Volpe 
2004, Serreau and Volpe 2004, McLaughlin 2004, Volpe 2005a), electroweak tests of 
the Standard Model (McLaughlin and Volpe 2004, Balantekin et al 2006a, Balantekin 
et al 2006b) as well as core-collapse Supernova physics (Volpe 2004, McLaughlin and 
Jachowitz 2005 and 2006). The physics case is also discussed in (Volpe 2005b, Volpe 
2005c, Volpe 2006). A crucial outcome of (Serreau and Volpe 2004) is that a small 
devoted storage ring appears as more appropriate for low energy applications. The 
feasibility study of this storage ring (concerning e.g. size, ion intensities, stacking 
method, space charge effects) is now ongoing (Chance and Payet 2005, Benedikt et 
al 2006). Here we summarize the status on the physics potential of low energy beta- 
beamgjj. 

3.1. Nuclear structure, neutrino-nucleus, nuclear astrophysics applications 

In (Volpe 2004) it has been pointed out that the availability of neutrino beams in the 
100 MeV energy range offers a unique opportunity to study spin-isospin and isospin 
nuclear excitations, for which nowadays little experimental information is available. 
Both the isospin and spin-isospin collective (and non-collective) modes are excited when 
a neutrino encounters a nucleus, due to the vector and axial- vector nature of the weak 
interaction. The (super-)allowed Fermi transitions - due to the vector current and 
therefore to the isospin operator - offer a well known example of such excitations at low 
momentum transfer. A precise knowledge of these nuclear transitions is essential for 
determining the unitarity of the CKM matrix, the analog of the MNSP matrix Eq.Q 
in the quark sector (see e.g. Hardy and Towner 2005). Another less known but still 
intriguing example is furnished by the allowed Gamow- Teller transitions - due to the 
axial-vector current and therefore to the spin-isospin operator - in mirror nuclei, which 

% Note that, after the proposals of pushing the ion Lorentz boosts to values higher than 60-100, the 
original scenario (Section II) is sometimes misleadingly referred to as "low energy" in the literature, to 
distinguish it from the medium and high-energy (Section IV). 
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Figure 16. Nuclear structure studies with low energy beta-beams: 
Contribution of isospin and spin-isospin nuclear states excited in the charged-current 
208 Pb(V e , e~) 208 Bi reaction (1CP 40 cm 2 ) for increasing neutrino energy, i.e. E Ue — 
15 MeV (up), 30 MeV (middle), 50 MeV (bottom). The histograms show the Isobaric 
Analogue State (J w = + ), the allowed Gamow- Teller (J 7 ' = 1+), the spin-dipole 
(J* = 0~,l - ,2-), as well as states of higher multipolarity (J 7r = 2+, 3 _ , 3+,4 - ,4+). 
For the latter no experimental information is available. Their contribution to the total 
cross section becomes significant when the impinging neutrino energy increases (Volpe 
2004). 



are a tool for the search of the possible existence of second-class currents in the weak 
interaction (Weinberg 1958, Wilkinson 2000). These terms transform in an opposite way 
under the 9- parity ^transformation as the usual vector and axial-vector terms, and 
are not present in the Standard Model. Because of their importance, both the allowed 
Fermi and Gamow- Teller transitions have been studied for a very long time in nuclear 
physics through beta-decay and charge-exchange reactions (for the Gamow- Teller ones). 
While the precision achieved in the description of the (super-) allowed Fermi transitions 
is really impressive, our understanding of the allowed Gamow- Teller transitions still 
requires the use of an effective axial-vector coupling form factor, to take into account 
the "quenching" of the predicted transitions, compared to the ones measured in beta- 
decay or charge-exchange reactions (see for example Osterfeld 1992). The still open 
"quenching" problem represents a limitation in our description of the weak spin-isospin 
nuclear response, in spite of the crucial role that it plays in various hot issues in nuclear 
astrophysics and in high energy physics. 

Little or no experimental information is available for the spin-isospin and isospin 
nuclear excitations such as the spin-dipole (J n = - ,l~,2~), or the states of higher 
multipolarity (e.g. J 71 = 2 + , 3~, 3 + , 4~, 4 + ). These states come into play when a 
weak probe transfers a finite momentum to a nucleus, like in muon capture or in 

+ The Q— parity transformation corresponds to the product of the charge-conjugation and of a rotation 
in isospin space. 
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neutrino-nucleus interactions. Supplementary information on the corresponding weak 
transition amplitudes can be furnished by electron scattering studies, which however 
explore nuclear excitations induced by the vector current only. As an illustrative 
example Fig{16] shows the contribution of spin-sospin and isospin transitions excited 
in neutrino scattering on lead, and their evolution when the neutrino energy increases. 
A quantitative estimate of the importance of such states is also gathered by computing 
the flux-averaged cross sections^, which are the relevant quantities for experiments. If 
one considers the neutrino fluxes corresponding to the decay-at-rest of muons, the spin- 
dipole states (J 7 " = 0~, 1~ , 2~) contribute by about 40 % in 12 C (Volpe et al 2000) and 
56 Fe (Kolbe and Langanke 1999), and by about 68 % in 208 Pb (Volpe et al 2002). The 
contribution from the states of higher multipolarity is about 5 % and 25 % in iron and 
lead respectively; while it amounts to about 30 % in carbon (Volpe et al 2000) and 
60 % in lead (Volpe et al 2002) if neutrinos are produced from pion decay-in-flight. 
Since low energy beta-beams have the specificity that the average neutrino energy can 
be increased by increasing the Lorentz boost of the ions, they constitute a promising 
tool for the study of these states, through a systematic study on various nuclear targets 
and different neutrino energies. Even though the measured cross sections are, in the 
majority of cases, inclusive, experimental information on these states can be extracted 
by changing the Lorentz ion boosts since different pieces of the nuclear response are 
important at different energies (Volpe 2004, Figure [TBI . 

A facility based on conventional sources (decay of pions and muons) offers another 
possible strategy to dispose of low energy neutrinos beams. The first feasibility studies 
show that the neutrino intensities are expected to be higher with conventional sources 
(Avignone F T et al 2000, Efremenko 2004) than with beta-beams. On the other hand, 
the neutrino beams obtained with the two options present complementary features, 
both for the flavour content and for the energy. Conventional sources provide us with 
neutrinos of different flavours and a Michel spectrum (from muon decay-at-rest) peaked 
at about 35 MeV and with 53.8 MeV maximum energy|]j. Beta-beams are pure in flavour 
and have the important advantage that the neutrino energy can be tuned by varying the 
acceleration of the ions. This is a particularly important feature for extracting nuclear 
structure information, as discussed above. In (McLaughlin 2004) neutrino scattering 
on lead is taken as an example to show another procedure to extract information 
on the different spin-isospin and isospin excitations, namely through a comparison of 
measurements with conventional beams and low energy beta-beams (Figure [T71) . In 
fact, the corresponding neutrino fluxes are in the same energy range for 7 = 7, but their 
shape and average energies are different. Besides, the measurement of the cross section 
without or with (one or two) neutrons can be used for the same purpose, in the specific 
case of the lead nucleus (McLaughlin 2004). 

* Flux-averaged cross sections are obtained by folding the cross sections with the neutrino flux of the 
source. 

jj See for example (Singh et al 2006b) for a recent calculation of neutrino-nucleus cross sections with 
various nuclear targets and conventional sources. 
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Figure 17. Total electron spectrum coming from charged current neutrino scattering 

on lead (solid line), with neutrinos coming from pion decay-at-rest (left) and from low 

energy beta-beams with 7 = 10 (right). The other curves present the contribution from 

the allowed Fermi and Gamow- Teller states (long dashed line), from the CP, 1~, 2 _ and 

the 2 + ,3~,3 + ,4~,4+ states (dotted line) (McLaughlin 2004). 

A better knowledge of the weak nuclear response can have a broad impact, 
since neutrino-nucleus interactions represent a topic of current great interest for 
neutrino physics and nuclear astrophysics. The motivations come from the need for 
a precise knowledge of neutrino detector responses, both in neutrino experiments 
and in core-collapse Supernovae observatories, and also from open issues in core- 
collapse Supernovae physics, such as unraveling the site for the nucleosynthesis of 
heavy elements (McLaughlin and Fuller 1995, Haxton et al 1997, Qian et al 1997, 
Borzov and Goriely 2000, Balantekin and Fuller 2003, Terasawa et al 2004). or from 
neutrino-nucleosynthesis (Woosley et al 1990, Heger et al 2005, Haxton 2004, Langanke 
and Pinedo 2004). However, the available experimental neutrino- nucleus scattering 
data in the relevant energy range is limited, since deuteron, carbon and iron are the 
only nuclei investigated so far. As a consequence one has to rely on the numerous 
theoretical predictions and on extrapolations to the nuclei and energies of interest. These 
calculations exploit a variety of approaches, among which Effective Field Theories, the 
Shell Model and the Random-Phase- Approximation, and the Elementary Particle Model 
(Kubodera and Nozawa 1994, Kolbe et al 2003). The cross sections estimates agree 
quite well at very low energies, where the nuclear response is dominated by the allowed 
Fermi and Gamow- Teller transitions. However, important discrepancies appear at higher 
neutrino energies, when the other nuclear excitations become important (Figures [T6l 
and HTj) . In this energy region, the calculations are in fact largely subject to nuclear 
structure uncertainties and model dependencies (e.g. treatment of the continuum, 
choice of the forces, higher-order correlations). The present discrepancies between the 
predicted and measured neutrino-carbon cross sections and between the calculations in 
the case of lead are talkative examples (Volpe 2005c) of the difficult theoretical task. 
Systematic neutrino-nucleus interaction studies performed with low energy beta-beams 
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Table 3. Neutrino-nucleus interaction rates at a low energy beta-beam 
facility (7 = 7— 14): The expected events per year (3.2 x 10 7 s) on several target 
nuclei obtained for 7 = 14 with Eqs. (|14M2"0")) are presented. The effect of the storage ring 
geometry is shown (S and D are the total and straight section lengths respectively). 
The detectors are located at 10 meters from the storage ring and have cylindrical 
shapes (i£=1.5 m and /i=4.5 m for deuteron, iron and lead, i?=4.5 m and h= 15 m 
for oxygen, where R is the radius and h is the depth of the detector). Their mass is 
indicated in the second column (Serreau and Volpe 2004). 



Reaction 


Mass 
(tons) 


Small Ring 

(5=450 m, D= 150 m) 


Large Ring 

(S=7 km, D=2.5 km) 


z/+D 


35 


2363 


180 


z/+D 


35 


25779 


1956 


z/+ 16 


952 


6054 


734 


Z7+ 16 


952 


82645 


9453 


z/+ 56 Fe 


250 


20768 


1611 


z/+ 208p b 


360 


103707 


7922 



offer a perfect tool to explore the nuclear response in this energy region in great detail, 
and put the theoretical predictions on a firm ground. 

The number of neutrino-nucleus scatterings expected at a low energy beta-beam 
facility has been calculated in (Serreau and Volpe 2004). The effects of the geometry of 
the storage ring has in particular been analyzed. Two sizes are considered: the one from 
the original scenario (Figure ED and a small storage ring, as the one envisaged for the 
future GSI facility (see GSI). The results of exact numerical calculations using Eqs.ffTH- 
l20l) are shown in Table [3]for several target nuclei. The detectors are those considered for 
the ORLAND facility (Avignone F T et al 2000). An analytical formula is derived to 
scale the rates in Table [3] to different storage ring sizes. Note that for a close detector, 
due to the flux anisotropy, the rates do not simply scale as S/D, as in the case of a 
far detector (like those considered in Sections II, IV and VI). It is clear that interesting 
neutrino-nucleus rates can be attained at such a facility. A crucial outcome of (Serreau 
and Volpe 2004) is that a small devoted storage ring is much more favorable for low 
energy applications. The physical reason is simple. The flux emittance is large because 
of the very low 7 values: only the ions which decay close to the detector contribute to 
the rates. The feasibility study of the small storage ring is now ongoing (Benedikt et al 
2006). 

A novel procedure to determine the response of a target nucleus in a supernova 
neutrino detector directly, through the use of low energy beta-beams, is pointed out in 
(Jachowicz and McLaughlin 2005 and 2006). It is shown that the cross sections folded 
with a supernova neutrino spectrum can be well reproduced by linear combinations of 
beta-beam spectra. This comparison offers a direct way to extract the main parameters 
of the supernova neutrino flux. The proposed procedure appears quite stable against 
uncertainties coming from the experiment, or the knowledge of the cross section, that 
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give rise to a "noise" in the expansion parameters. 

Finally, it has been recently pointed out (Volpe 2005a) that neutrino-nucleus 
interactions are also important for the search of neutrinoless double-beta decay in nuclei. 
In fact, by rewriting the neutrino exchange potential in momentum space and by using 
a multiple decomposition, the two-body transition operators, involved in the former, 
can be rewritten as a product of the one-body operators involved in neutrino-nucleus 
interactions (except for the short range correlations as well as possible phases present in 
the two-body process). Neutrino-nucleus scattering data offer a potential new constraint 
for the predictions on the neutrinoless double-beta decay half-lives. At present these 
calculations suffer of important discrepancies for the same candidate nucleus. Beta- 
decay (Muto et al 1989, Aunola and Suhonen 1996), muon capture (Kortelainen and 
Suhonen 2002 and 2004), charge-exchange reactions (Akimune 1997, Bernabeu et al 
1988) and double-beta decay with the emission of two neutrinos (Rodin et al 2003) 
have been used to constrain the calculations so far. Neutrino-nucleus measurements 
would have the advantage that, if both neutrinos and anti-neutrinos are available, the 
nuclear matrix elements involved in the two branches of neutrinoless double-beta decay 
- from the initial and the final nucleus to the intermediate one - can be explored. 

3.2. Fundamental interaction studies 

Three applications for fundamental interaction studies of low energy beta-beams have 
been discussed so far: the measurement of the Weinberg angle at low momentum transfer 
(Balantekin et al 2006a), a CVC test with neutrino beams (Balantekin et al 2006b) and 
a measurement of the neutrino magnetic moment (McLaughlin and Volpe 2004). The 
measurement of the Weinberg angle represent an important test of the electroweak 
theory. Several experiments at different Q 2 exist, namely the atomic parity violation 
(Bennett and Wieman 1999) and Moller scattering at Q 2 = 0.026 GeV 2 (Anthony et 
al 2005) which combined with the measurements of sin 2 9w at the Z° pole (ALEPH et 
al 2005), are consistent with the expected running of the weak mixing angle. However, 
recent measurement of the neutral- to charged-current ratio in muon anti-neutrino- 
nucleon scattering at the NuTEV experiment disagrees with these results by about 3 a 
(Zeller et al 2002). A number of ideas were put forward to explain the so-called NuTEV 
anomaly (Davidson et al 2002, Loinaz et al 2004, Miller and Thomas 2005, Giunti and 
Laveder 2002). However, a complete understanding of the physics behind the NuTEV 
anomaly is still lacking; probing the Weinberg angle through additional experiments 
with different systematic errors would be very useful. The possibility of using a low 
energy beta-beam facility to carry out such a test with neutrino-electron scattering at 
low momentum transfer, i.e. Q 2 = 10~ 4 GeV 2 is investigated in (Balantekin et al 2006a). 
In particular it is shown that since the neutrino flux and average energy are well known 
in the case of beta-beams, the number of counts is in principle sufficient to extract 
information on the Weinberg angle. A fully efficient 1 kton Cerenkov detector is located 
at 10 meters from a small storage ring having 1885 m total length and 678 m straight 
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Figure 18. A Weinberg angle measurement with low energy beta- 
beams: A% 2 with 7 = 12 (dotted line), with 7 = 7,12 (broken line), and with 
7 = 7, 8, 9, 10, 11, 12 (dash-dotted line). The results were obtained considering a one 
year (3 x 10 7 s) measurement duration at each 7, and for a helium-6 intensity at 
the storage ring of 2.7 x 10 12 ions/s (Chance and Payet 2005). The count number 
error was considered to be purely statistical. The la (Ay 2 = 1) relative uncertainty 
in the Weinberg angle is 15.2% for 7 = 12, 12.3% for 7 = 7,12, and 7.1% for 
7 = 7, 8, 9, 10, 11, 12 (Balantekin et al 2006a). 



sections, in which the expected intensities for 7 = 7 — 14 are 0.5X10 11 helium-6/s and 
2.7xl0 12 neon-18/s (while the intensities at production are the same as in Tabled]). 
These numbers are the outcome of a preliminary feasibility study of the small storage 
ring (Chance and Payet 2005). The (anti)neutrino on electron events are identified by an 
angular cut. The background from neutrino-proton scattering is suppressed by the use of 
Gadolinium (Beacom and Vagins 2004). Figure [181 shows the improved precision one gets 
by combining measurements with different gamma values. The results corresponds to 
one (3xl0 7 s) year measurement at each gamma and a purely statistical total count error. 
Figure [19] shows the precision with which the Weinberg angle can be measured after the 
inclusion of both statistical and systematic errors. In particular, if the systematic error 
can be kept below 10 %, a precision of 10 % seems to be within reach at a beta-beam 
facility. 

The CVC hypothesis connects weak and electromagnetic hadronic currents. Several 
tests of CVC have been performed in the past (see e.g. Thomas 1996), concerning, in 
particular, the vector form factor, through super-allowed nuclear beta-decays studies 
(see e.g. Hardy and Towner 2004). Verifying that the CVC hypothesis correctly 
predicts tensor terms - often referred to as weak magnetism - is of fundamental 
importance (Deutsch et al 1977). So far, this contribution to the weak currents has 
been tested in an experiment involving the beta-decay of Gamow- Teller transitions in 
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Figure 19. One sigma uncertainty in the Weinberg angle as a function of the 
systematic error at each 7 for 7 = 12 (dotted line), for 7 = 7, 12 (broken line), and for 
7 = 7, 8, 9, 10, 11, 12 (dash-dotted line). The helium-6 intensity at the storage ring is 
2.7 x 10 12 ions/s on the left panel (Chance and Payet 2005), and with an increased 
intensity of 8.1 x 10 12 ions/s on the right panel. In both cases, the measurement 
duration at each 7 is one year (3 x 10 7 s) (Balantekin et al 2006a). 
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Figure 20. A CVC test: The v e + p — > e+ + n reaction is used to test the 
weak magnetism form factor using neutrinos from low energy beta-beams and a water 
Cerenkov detector. The figure shows the ratio of the first {A-\) and second (A 2 ) order 
Legendre polynomial coefficients to the zeroth order one (Aq, the number of events) 
as a function of the ion Lorentz boost. These expansion coefficients determine the 
angular distribution of the emitted positrons with (solid is for A\/Aq\ dash-dotted is 
for A2/A0) and without (dashed is for A\/Aq; dotted is for A2/A0) the weak magnetism 
term. The overall effect of the weak magnetism term on the angular distribution is 
considerably larger than on the total number of events and increases with the Lorentz 
ion boost. Note that when the weak magnetism is considered, its value is set to its 
CVC value (Balantekin et al 2006b). 



mirror nuclei in the A=12 triad (Lee et al 1963, Wu 1964, Lee and Wu 1965). A test 
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Figure 21. Same as for the previous Figure, but for A% 2 obtained from the angular 
distribution of events (see text) in the cases when the total error is purely statistical 
(solid), and when one includes 2% (dashed), 5% (dash-dotted) and 10% (dotted) 
systematic errors. The la (A^ 2 = 1) relative uncertainty in // p — //„ is 4.7%, 5.6%, 
9.0% and more than 20%, respectively. These results were obtained considering 7=12 
(Balantekin et al 2006b). 



based on neutrino-nucleon collisions at low momentum using low energy beta-beams is 
proposed in (Balantekin et al 2006b). This would have, in particular, the advantage 
that there is no uncertainty coming from nuclear structure calculations. It is important 
to keep in mind that, in addition to providing important tests of the Standard Model, 
neutrino-nucleon reactions play a crucial role in understanding the dynamics of the 
core-collapse supernovae (Balantekin and Fuller 2003, Horowitz 2001), the yields of 
the r-process nucleosynthesis that could take place in such environments (Meyer et al 
1998) and also contribute to the energy transfer (from the accretion- disk neutrinos to 
the nucleons) in gamma-ray bursts models (Ruffert et al 1997, Kneller et al 2004). 
Finally, understanding the subtleties of the neutrino-nucleon interactions are crucial to 
the terrestrial observation of neutrino signals (Vogel and Beacom 1999, Beacom et al 
2002). In (Balantekin et al 2006b) the sensitivity to the weak magnetism term that can 
be achieved in the v e + p — > e + + n reaction, both with the total number of events and 
with the angular distribution of the emitted positrons, is studied. For this purpose the 
same setup, as for the measurement of the Weinberg angle is taken, and the interaction 
of anti-neutrinos on protons in a water Cerenkov detector is considered. The angular 
distribution can be expanded in Legendre polynomials whose coefficients are shown in 
Figure |20| as a function of the Lorentz ion boost. While the zeroth order coefficient is 
nothing but the total number of events, the first order measures the strength of a linear 
term which favors forward peaked events. The results show that, when systematic 
errors are taken into account, the angular distribution is a much better tool than the 
total number of events to extract information on the weak magnetism form factor. In 
particular, if those errors are kept below 5 %, a one year measurement of the weak 
magnetism is possible at a 1 a level of 9 %, if the ions in the storage ring are boosted 
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to 7 = 12. An even better measurement is expected if the ions are boosted to Ys larger 
than 12, because of the increasing importance of the weak magnetism contribution with 
the impinging neutrino energy. This way of probing the weak magnetism form factor 
at low momentum transfer constitutes a new test of the Conserved Vector Current 
Hypothesis. 

It is interesting to note that one might use the ion decay at rest as an intense 
neutrino source, in order to explore neutrino properties that are still poorly known, such 
as the neutrino magnetic moment (McLaughlin and Volpe 2004). Direct measurements 
to achieve improved limits are precious, since the observation of a large magnetic moment 
points to physics beyond the Standard Model. Note that such an application does not 
require any acceleration of the ions. Once produced, the ions are fired to a target inside 
a Air detector. The measurement of the (anti)neutrino interaction with the electrons, 
as a function of the electron recoil, is then used to set limits on the neutrino magnetic 
moment. Current bounds come from direct measurements using reactor neutrinos, and 
are in the range \x v < 0.9 — 4 x 10 -10 /i£ (Daraktchieva et al 2005, Li et al 2003, Reines 
et al 1976, Vogel and Engel 1989, Vidyakin et al 1992, Derbin et al 1993). From 
solar neutrino-electron scattering at Super-Kamiokande a limit of < 1.5 x 10~ 10 /i£ at 
90 % C.L. has been obtained (Beacom and Vogel 1999). Upper limits in the range 
10 -11 — 10~ 12 [iB are also inferred, from astrophysical and cosmological considerations 
(Raffelt 1996), the exact values being model dependent. The prospects from using low 
energy beta-beams to improve the direct bounds has been studied in comparison with 
reactor neutrinos and a very intense tritium source (McLaughlin and Volpe 2004). While 
the advantage of using beta-beams is that the neutrino flux is perfectly known, the main 
limitation for this application is clearly the intensity of the ions. Figure [22] shows the 
electron-neutrino scattering events in the range of 0.1 MeV to 1 MeV and 1 keV to 10 
keV respectively. The results shown are obtained for an (increased) intensity of 10 15 u/s 
(i.e. 10 15 ions/s). If there is no magnetic moment, this intensity produces about 170 
events in the 0.1 MeV to 1 MeV range per year and 3 events in the 1 keV to 10 keV 
range per year. These numbers increase to 210 and 55 respectively in the case of a 
magnetic moment of 5 x lO -11 /^- In conclusion, by measuring electron recoils in the 
keV range with a beta-beam source one can, with a sufficiently strong source, have a 
very clear signature for a neutrino magnetic moment of few lO -11 ^. A detailed study 
of the background, at these small electron recoil energies, and a precise knowledge of 
the ion-induced background, are now needed. 

4. Medium and high energy scenarios 

The idea of pushing the Lorentz ion boosts to higher values than the ones proposed in 
the original scenario (Zucchelli 2002) is first proposed in (Burguet-Castell et al 2004). 
In particular, the sensitivity to CP violation of the original scenario (Section II) is 
compared to two new baselines where the ions are accelerated at either 7 = 350/580 
for helium-6/neon-18, and sent to a distance of L = 720 km, or at 7 = 1500/2500 




Figure 22. MEASUREMENT OF THE NEUTRINO MAGNETIC MOMENT WITH 

LOW energy beta-beams: Number of neutrino-electron scattering events 
from helium-6 ions, produced at the rate 10 15 per second, and sent inside a 
4tt detector of 10 m in radius. Similar results are obtained with neon-18 ions. 
(Note that for this application ion acceleration is not needed.) The symbols 
show the number of scatterings if the neutrino has a magnetic moment of 
fj,„ = W~ 10 fiB (diamonds), fi v = 5 x 10~ 11 /j lB (stars) and [i v = 10 -11 /i£ 
(triangles). The histogram shows the results for a vanishing neutrino magnetic 
moment (McLaughlin and Volpe 2004). 



and sent to L = 3000 km. The former case would require a refurbished SPS to reach 
such accelerations at CERN, while the water Cerenkov detector detector could, for 
example, be located at the Gran Sasso Laboratory. The latter case would necessary 
need something like LHC and a baseline such as the CERN-Canary islands with e.g. a 
40 kton tracking calorimeter. If the baseline is adjusted appropriately one can remain 
close to the atmospheric oscillation maximum. Working at higher boosts has essentially 
two advantages. The appearance oscillation signal increases linearly with gamma (and 
then saturates for 7 « 400 in a water detector), which results in a gain in sensitivity. 
Besides, some ambiguities on the neutrino parameters due to the degeneracies can be 
disentangled, both by using the neutrino energy reconstruction in the detector, and by 
exploiting matter effects (that become significant for baselines of 0(1000) km, Figure 
[5]). Ref.(Burguet-Castell et al 2004) has triggered a number of studies (Terranova et al 

2004, Burguet-Castell et al 2005b, Huber et al 2006, Agarwalla et al 2005). The medium 
and high energy options have also been discussed in (Terranova 2005, Gomez-Cadenas 

2005, Migliozzi 2005). In (Donini et al 2005a and 2006) the medium 7 (350-580) case 
is analyzed exploiting a CERN-Gran Sasso Laboratory baseline with a non-magnetized 
iron detector. It is shown that the CP violating phase can be determined down to 30° 
(99 % C.L.), if #13 > 3°. In (Agarwalla et al 2005) the high energy option is discussed 
with a baseline from CERN to Rammam or Singara where an iron calorimeter detector 
could be located. Performances can be found in the literature for 7 = 150 (Burguet- 
Castell et al 2005), 200 (Huber et al 2006), 350 (Burguet-Castell et al 2005), 500 (Huber 
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Figure 23. Sensitivity curves (99 % C.L.) for 6 (upper) and #13 (lower), as a function 
of the ion boost. The CERN-Frejus baseline with a 440 kt water Cerenkov detector are 
considered. The (anti)neutrino emitters run at the same gamma. The stars present 
the results obtained for the original option where the ions run at 7 = 60, 100 (Section 
II) (Burguet-Castell et al 2005). 



et al 2006, Terranova et al 2004), 1000 (Ruber et al 2005a), 2000 (Burguet-Castell et al 
2004) and 2488 (Terranova et al 2004). Precision tests of several observables such as the 
weak mixing angle, the neutrino magnetic moment and potential leptonic Z' couplings, 
exploiting neutrino-electron scattering and the medium energy (7 = 500) beta-beam 
are discussed in (De Gouvea and Jenkins 2006), in comparison with nuclear reactors, 
neutrino factories and conventional beams. 

In (Burguet-Castell et al 2005) a detailed analysis of the sensitivity with a Lorentz 
ion boost of 7 = 150 is made. This choice corresponds to the maximum acceleration 
that can be achieved using the existing CERN-SPS (Lindroos 2005). In particular, the 
following three issues are addressed: i) a sensitivity study of the CERN-Frejus baseline 
with a megaton water Cerenkov detector and 7 > 100; ii) the optimal baseline at this 
7 and with such setup; in) the (dis)advantages of using different values of 7 for the 
(anti) neutrino emitters. The intensities used are the same as in the original scenario, 
namely 2.9xl0 18 6 He and l.lxlO 18 18 Ne, while the measurement period is 10 years. 
Note that the water Cerenkov detector considered is suitable for quasi-elastic events. 
Therefore the sensitivity improves until the inelastic cross section starts to dominate 
(for 7 > 400). The angular and energy cuts as well as the time structure of the ion 
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99% CL sensitivity to CP violation, for 4.4 Mton y 99% CL sensitivity to 6 13 , for 4.4 Won y 




sin 2 (28, 3 ) 5{deg) 

Figure 24. Exclusion plots (99 % CL.) for 8 (left) and 6*13 (right) for the original 
scenario in comparison with the options 7 = 120, 150,350 and L = 130,300, 730 km 
for the source-detector distance, respectively. The solar and atmospheric oscillation 
parameters are fixed to the present best values (Burguet-Castell et al 2005). 




Figure 25. The figure shows a comparison of the CP violation sensitivity for the 
normal (left) and inverted (right) mass hierarchy as a function of the true values for 
sin 2 and Sep at the 3 a CL. A comparison is made between beta-beams with 
different ion boosts and detector technologies (WC and TASD) on one hand, and 
neutrino factories (NF@3000 km) and super-beams (T2HK) on the other hand. On 
the y-axis is shown the values of 5c p "stacked" to the CP fraction (Huber et al 2006). 



bunches help keeping the atmospheric background at a very low level, while the beam- 
induced neutral current single-pion production, that can mimic the appearance signal, 
is manageable even for 7 = 100. 

Figure [23] shows the sensitivity on 5 and #13 as a function of the ion boost for the 
CERN-Frejus baseline. It is clearly seen that by using 7 = 100 the increased event 
rates imply a definite gain in the 6 sensitivity compared to the lower gamma values 
first envisaged, while its behavior keeps rather flat in the range 100-150 (since for these 
cases the neutrino flux at lower energies contributing to the event rate is practically the 
same). The study of the optimal source-detector L distance, for the maximum SPS ion 
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Comparison of discovery reaches (3cr) 
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Figure 26. A comparison of the discovery reach at the 3cr confidence level. The bars 
represent the best case (left end), "typical case" (middle line), and worst case (right 
end) in Sep, where the respective sensitivities are computed with the CP fractions in 
the plot legend. Note that for CP violation, a CP fraction of one can never be achieved 
since values close to or it cannot be distinguished from CP conservation. Within each 
category, the most competitive setups from this study are compared. For all shown 
sensitivities, a normal mass hierarchy is assumed (Huber et al 2006). 



acceleration of 150, points to L — 300 km, i.e. the distance corresponding to the first 
atmospheric neutrino oscillation maximum. 

A comparison of the 5 and #13 exclusion plots for the original scenario, the 
7 = 120, 150 with L = 130, 300 km, and 350 with and 730 km from (Burguet-Castell 
et al 2004), is presented in Figure EH Note that the energy dependence of the signal 
can be well reconstructed since the neutrino energies are above the Fermi momentum 
of the target nuclei (this is not the case for the original scenario). This information can 
be used at 7 = 150 (350) and L = 300 (700) to resolve the intrinsic degeneracies. A 
major issue for such a baseline is to find a location to house the detector, the Canfranc 
Underground Laboratory being now discussed as a possible option. 

An interesting and detailed analysis of the CV violation sensitivity is made in 
(Huber et al 2006), for Lorentz ion boosts values ranging from 100 to 3000, and is also 
discussed in (Winter 2006). Note that here the performance at high gammas is studied 
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considering that the intensities are reduced (and not constant as in previous studies) 
for very high gammas. While a 500 kt water Cerenkov detector is considered for a 
boost of 200, a 50 kt totally active scintillator detector (TASD) is taken for higher ion 
accelerations. The whole study is performed by using the GLoBES software (Huber et al 
2005a) which allows, in particular, to include degeneracies and correlations (see Section 
V) . A detailed comparison with neutrino factories with a 50 kt magnetized iron detector 
and very intense super-beams (as those that will be available at T2HK) is made (Figure 
l25i) . Figure [261 presents the main results of this work as far as the third neutrino mixing 
angle, the Dirac phase and the normal mass hierarchy are concerned. One can see 
the excellent performance of the 7 = 200 setup compared both to super-beams and to 
neutrino factories. Such a discovery reach is achieved by including the second oscillation 
maximum to disentangle correlations and degeneracies. This also depends on the choice 
of approximately equal neutrino and anti-neutrino event rates for all setups. On the 
other hand the higher gammas appear as competitive options to a neutrino factory. In 
conclusion, the discovery reach for 7 > 150 is very promising as well. However, at least 
for the moment, it remains more speculative, since it is neither based on an existing 
accelerator complex nor on a robust estimation of the ion intensities. These scenarios 
clearly deserve further investigations. 

5. The issue of the degeneracies 

A characteristic feature in the analysis of future neutrino accelerator experiments is the 
presence of the so-called "parameter degeneracies", which is the appearance of various 
disconnected regions in the multi-dimensional oscillation parameter space, besides those 
corresponding to the true solutions. The appearance of such degeneracies is due to 
both the inherent three flavour structure of the oscillation probabilities, and some still 
unknown neutrino oscillation parameters. The presence of the latter degeneracies in the 
analysis of a beta-beam experiment depends on our knowledge of neutrino oscillation 
properties at the time such measurement is performed. These degeneracies are of four 
different kinds, namely the intrinsic, the sign, the octant and the mixed degeneracies 
(or clones). The classification corresponds to the following origins: 

• for an appearance experiment the oscillation probability Eq.( fT3l) obtained for 
neutrinos at fixed baseline and energy has a continuous number of solutions 
(S, 613). If the same experiment is performed for both neutrinos and anti-neutrinos 
with the same energy, only two disconnected solutions are left that satisfy the 
equation P a a(5, 613) = P^((5, ^13) (where ± refers to neutrino and anti-neutrino 
respectively). These two regions correspond to the true solution and the intrinsic 
clone (Burguet-Castell et al 2001). 

• The sign of Am^ is still unknown. The positive and negative solutions give two 
different regions called the sign clones (Minakata and Nunokawa 2001). 

• The two octants 823 > it /4 and 623 < vr/4 are possible, giving in general two possible 
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Figure 27. The figure shows the 90% CL contours for 6>i 3 = 1°, 7° and 5 = 90° (light 
red circle), after a 10 year /3-Beam (upper), a 2+8 year Super-Beam (middle) and a 
combined measurement (lower figure). The different clone regions are shown (see text): 
the intrinsic (solid), the sign (dashed), the octant (dot-dashed) and the mixed (dotted 
line) degeneracies. The theoretical clone locations are presented as well (dark black 
dots) (Donini et al 2005b). 
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Figure 28. CP discovery potential in the (613, S) plane at 99 % C.L. for a super-beam 
(left) and a beta-beam (right), considering the appearance channel only (dashed line), 
and a combination of the appearance and the disappearance channels (solid line) at the 
same facility. A 2% systematic error in the disappearance channel and a 5% systematic 
error in the appearance channel have been considered for both facilities (Donini et al 
2005c). 

solutions for S and #13. This is called the octant degeneracy (Fogli and Lisi 1996). 

• A "wrong" sign and a "wrong" octant give rise to a new degeneracy called the 
mixed clone. 

The ensemble of these four regions for the true solution and the intrinsic clone gives eight 
different solutions also called the eightfold degeneracy (Barger et al 2002). In (Donini 
et al 2004) it is pointed out that the clone solutions can be significantly different if 
determined from the oscillation probability, or from the number of events, showing that 
the degenerate solutions should be computed from the number of events, which is the 
physically measured quantity. Analytical formulae to determine the clone locations 
as a function of the real solutions are given. The effect of possible combination of 
super-beams and beta-beams (Section II) as well as of a neutrino factory with super- 
beams and/or beta-beams is discussed. A detailed analysis of the eightfold degeneracy 
is performed for the first time in (Donini et al 2005b). It is shown that combining 
a super-beam and a beta-beam does not solve the problem of degeneracies, since, at 
the neutrino energies involved, both are counting experiments, tuned at the same L/E 
value. Some clone regions are eventually reduced (Figure [27]). Different strategies can 
be followed to resolve the degeneracy issue, such as the use of spectral analysis (Burguet- 
Castell et al 2001), combination of various experiments (Minakata et al 2003), and/of 
different channels (Donini et al 2002). In (Donini et al 2005c) the combination of the 
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Figure 29. Impact of the degeneracies on the CP discovery potential for the 7 = 100 
beta-beam. We show the sensitivity at 3a (A% 2 > 9) computed for 4 different 
options about the true parameter values: Normal: sign(Am^ 1 ) — 1, #23 = 40°; 
Octant: sign(Am§ 1 ) = 1, 8 23 = 60°; Sign: s^7i(Am§ 1 ) = -1, 9 23 = 40°; Mixed: 
sz(7n(Am| 1 ) = —1, #23 = 60°. Dotted curves are computed neglecting degeneracies 
(Campagne et al 2006). 

appearance and disappearance channels is explored both for the original beta-beam 
and for a super-beam, showing that while there is no improvement for the former as 
soon as some level of systematic error is added, degeneracies are reduced for the latter. 
Figure [28] presents the discovery potential for each of these options taken alone, taking 
fully into account the degeneracies. Note that the CP sensitivity contours have only an 
approximate symmetry for \b\ > 7r/2 and \5\ < tt/2, while there is no symmetry between 
positive an negative values of 5. The impact of the degeneracies on the CP discovery 
reach for the 7 = 100 option is also shown in Figure [29] (Campagne et al 2006). A 
discussion on the issue of degeneracies can be found in (Rigolin 2005 and 2006) as well. 
In (Donini and Fernandez-Martinez 2006), and (Fernandez-Martinez 2006), the idea of 
alternating the ions, i.e. 6 He and 18 Ne with 8 B and 8 Li, boosted at the same 7 of 100, 
is proposed to reduce the degeneracies. Indeed, high 8 B and 8 Li production rates could 
be possible thanks to a new production method called "ionsation cooling" (Rubbia et 
al 2006). Since the Q- values for the two sets of ions are very different, the neutrino 
fluxes would be peaked at the second oscillation maximum, for the former, and at the 
first oscillation maximum, for the latter, at a fixed baseline of 700 km. As a result, it is 
shown that such a measurement, even though less sensitive to #13, reduces the eightfold 
degeneracy to a two- fold one, thanks to the different L/E dependence. This option is 
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Figure 30. Impact on the degeneracies of a combined accelerator and atmospheric 
data analysis, in the case of a 7 = 100 beta-beam and a super-beam (SPL) at CERN 
on one hand, and of T2K phase II to Hyper-K (Itow et al 2001) on the other hand. The 
figure shows the allowed regions in sin 2 20i3 and (5cp for accelerator data alone (contour 
lines) and accelerator plus atmospheric data combined (colored regions). Solution for 
the true/ wrong mass hierarchy (octant of 23 ) are indicated with H tr / wr (O tr / wr ). The 
true parameter values are i5cp = — 0.85-7T, sin 2 26*13 = 0.03, sin 2 6*23 = 0.6 (Campagne 
et al 2006). 



particularly well suited if 6*13 is measurable at T2K-I. 

Another very interesting idea proposed recently to resolve the open problem of 
the degeneracies is to combine accelerator experiments with data from atmospheric 
neutrinos provided by the same detector (Huber et al 2005b). In fact atmospheric 
neutrinos are sensitive to #13 and the mass hierarchy due to earth matter effects in the 
e-like events. In addition effects from the solar parameters and Am^ on e-like events 
in the sub-GeV energy range provide a sensitivity on #23 . Due to these effects the sign 
and octant degeneracies can be broken. A first application is done in (Campagne et al 
2006) which shows that the true solutions are identified by the combined accelerator 
and atmospheric data both for a beta-beam and a super-beam (SPL) at CERN (Figure 
1301) . In particular, the hierarchy can be identified at 2a C.L. from a combined analysis 
if sin 2 9 13 > 0.03 — 0.05 for the beta-beam and the super-beam while the octant can 
be determined through the super-beam and the atmospheric data. Let us recall that 
these experiments alone have no sensitivity at all on the octant and sign issues. This 
important synergy between accelerator and atmospheric data is clearly very promising 
for the identification of the true value of the third mixing angle as well as of the CV 
violating phase in a beta-beam experiment. 

6. Monochromatic neutrino beta-beams 

An exciting variant to the original beta-beam idea has been proposed, where 
monochromatic neutrino beams are produced through the electron capture of radioactive 
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Table 4. Decay properties of some rare-earth nuclei. The half-life and the branching 
ratios for electron-capture relative to beta-decay are given in the first columns, while 
the last column presents the neutrino energy in the rest frame, defined as the Q^-value 
minus the energy at which the Gamow- Teller resonance is located (see text). The given 
ions are example candidates for a beta-beam facility producing monochromatic beams 
through electron capture (Bernabeu et al 2005a). 



Decay 


t 1/2 EC/ 13+ (%) £ (keV) 


148 Dy ^148 Xb 
150 Dy ^150 Tb 

152 Tm 2- ^ 152 Er 
i50 Ho 2 - ^i50 Dy 


3.1 m 96/4 2062 

7.2 m 99.9/0.1 1397 
8.0 s 45/55 4400 
72 s 77/33 3000 



ion beams (Bernabeu et al 2005a, Sato 2005). The baseline of such a facility would be 
the same as for the original beta-beam facility (Zucchelli 2002) (Fig. [8]). The production 
of such beams depends on the existence of ions decaying through electron capture with a 
short enough lifetime for electron capture to occur. In fact one needs to keep one electron 
bounded to the ion and therefore to produce, accelerate and store partly charged ions. 
Since, even in a very good vacuum conditions, a few remaining atoms make the ions 
quickly lose the electron, it is crucial that their lifetime is short enough to have high 
neutrino intensities. Good candidates are offered by some recently discovered nuclei 
far from the stability line, which have a fast decay through an allowed transition to a 
Gamow- Teller resonance in the energetically accessible region. Examples of such nuclei 
are given in Table HI The monochromatic neutrino flux produced by electron capture 
and arriving in a detector at a given distance L is 

i^^^L HE - 27Eo) (22) 

with 7 >> 1, Ni ons is the total number of decaying ions, T v /T is the branching ratio for 
electron capture (Table HJ). The quantities E and E = E /[y(l — (3 cos 6)] represent the 
energy in the rest and in the laboratory frames respectively, with the angle 9 being the 
deviation of the neutrino with respect to the prolongation of the long straight section 
of the decay ring. 

There is a major difference in the use of neutrino beams coming from electron 
capture and from beta-beams in the search for CV, since only electron neutrino beams 
can be produced in the first case. As a consequence, the study of CV violation cannot 
be made through a comparison of neutrino versus anti-neutrino oscillations Eq. (fT0l) . 
but by performing oscillation measurements at different ion boosts. In fact, for a fixed 
distance between the source and the detector, the appearance probability for u e — > 
being energy dependent, the effect of the 5 phase is to shift the oscillation interference 
pattern, while the amplitude of the oscillation probability is determined by the mixing 
angle #13 (Figure EI]) (Bernabeu et al 2005a). As a consequence the values of these 
unknown parameters can be extracted by comparing the v e — > oscillation probability 
at different neutrino energies (ion boosts). The main advantage of the electron capture 
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Figure 31. Appearance probability for neutrino oscillations in the three flavour 
formalism and a non-zero Dirac 5 phase, as a function of neutrino energy, with a 
distance between the source and the detector of 130 km. The three curves correspond 
to three different values of the phase S, while the vertical lines show the two ion Lorentz 
boosts used in the simulation (Bernabeu et al 2005a) . 



case over beta-beams is that in the former case the energy of all neutrinos can be peaked 
at the energy (ies) of interest, to match the the oscillation maxima; while in the latter 
case, because of the broad spectrum (Figj6]), for a given intensity, many neutrinos will 
be in an energy range which is less sensitive to the phase and/or where the cross section 
is low. 

The sensitivity to CP violation has been studied in (Bernabeu et al 2005a), also in 
combination with a f3~— beam (Bernabeu et al 2006). Figure [321 shows the sensitivity in 
the (5, #13) plane for four chosen values of these parameters. The source is of 10 18 150 Dy 
ions/year, for a running time period of 10 years. The ions run 5 years at the Lorentz 
boost 7 = 195 - the maximum achievable at the CERN's SPS (FigjHJ) - and 5 years at 
7 = 90, chosen to be as different as possible from the first one and high enough to avoid 
backgrounds. The water Cerenkov detector has 440 kton fiducial mass and is located 
at a distance of 130 km (the CERN-Frejus distance). 

The possibility of using monochromatic electron neutrino beams to perform CP 
violation searches appears very promising. A detailed study concerning the ion 
production rates, the target and ion source design, the accumulation scheme as well 
as possible vacuum improvements is now required, to definitely state the achievable 
flux. 
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Figure 32. Sensitivity to the CV violating phase S and the still unknown mixing 
angle # 13 . The results are obtained considering neutrino beams from electron-capture 
of partly charged 150 Dy ions, boosted at 7 = 90, 195. The different curves show the 
confidence levels for the assumed values (8°, 0), (5°, 90°), (2°, 0) and (l°,-90°). The 
details on the used setup can be found in the text (Bernabeu et al 2005a). 



7. Conclusion and perspectives 

The recent discoveries in the field of neutrino physics open important new developments. 
Some require in particular pure, intense and well controlled neutrino beams, such as 
those furnished by "beta-beams". This novel method exploits the intense radioactive 
ion beams under study, and offers important synergies with the nuclear physics 
community. A new facility based on beta-beams has been proposed with the main 
goal of addressing the crucial issue of the existence of CV violation in the lepton 
sector. In the original scenario information on the CV violating phase is extracted 
from a comparison of neutrino versus anti-neutrino oscillations. The neutrino beams 
are produced at CERN, and sent to a gigantic Cerenkov detector located in the Frejus 
Underground Laboratory, 130 km from CERN. The first steps of the facility are very 
similar to the EURISOL project, while the acceleration to GeV energies uses already 
existing accelerator infrastructures (the PS and SPS) which implies important savings 
in the cost. The ions are then stored in a large storage ring that needs to be built. 
This scenario has also the significant advantage that it uses reasonable extrapolation of 
existing technology. 

While a beta-beam alone already has a very interesting discovery potential, it can 
be pushed even further if a very intense conventional beam is fired to the same detector 
as well. This combination also offers the possibility to perform studies of the T and 
CVT symmetries. On the other hand, it has recently been pointed out that if the CVT 
symmetry is assumed, the neutrino measurements with a beta-beam and a super-beam 
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in the T channel can be used to reduce the measurement time from 10 to 5 years. The 
achieved sensitivity is very close to that of a single experiment running for 10 years in 
both the neutrino and the anti-neutrino modes. 

The discovery reach of the original scenario is being explored in great detail. Various 
options are considered, such as running the neutrino and antineutrino emitters either 
together or independently in the storage ring, and/or varying the Lorentz ion boost 
from 60 to the maximum acceleration achievable in the SPS. (If at 7 of about 60 the 
experiment is a counting experiment only, already at about 100 some information on 
the energy of the events can be exploited.) This facility is well positioned to explore 
values between 10~ 3 to 10~ 4 for the third still unknown neutrino mixing angle #13. For 
the 7 = 100 option the maximal violation corresponding to 5 = n/2 and 3ir/2 can be 
discovered down to sin 2 2^13 = 6(8) 10~ 4 (99 % CL.) for the beta-beam (super-beam); 
while the best sensitivity is obtained for sin 2 29 13 > 10~ 2 since at this value CV violation 
can be established for 73 % (75 %) for all 5 values with a beta-beam (super-beam). The 
comparison with e.g. T2K phase II to Hyper-K in Japan - a project of similar size and 
timescale - shows a very similar discovery reach. Besides, the MEMPHYS detector has 
the same physics potential for proton decay and for (extra-) galactic supernova neutrinos 
as UNO in the US and Hyper-K. 

A variant of the original beta-beam idea consists in using ions decaying through 
electron-capture instead of beta-decay. This option has the advantage that the neutrino 
beams are mono-chromatic. (No anti-neutrino beams can be produced.) Information 
on the S phase is extracted in this case by performing experiments at different ion 
boosts, since the phase introduces shifts in the appearance oscillation patterns. The 
recent discovery of exotic nuclei that have fast decay through a strong Gamow- Teller 
transitions offer good potential candidates for the neutrino emitters. The required 
baseline is here the same as in the original beta-beam scenario. A specific feasibility 
study is now required as far as the achievable intensities and vacuum issues are concerned 
for example, since not fully stripped ions need to be accelerated and stored. Also, the 
sensitivity studies are to be pushed further. 

Medium and high energy scenarios have been proposed as well, where the ion 
boosts are of several hundreds and thousands respectively. Such energies require new 
baseline scenarios: new accelerator infrastructures, such as a refurbished SPS (for the 
medium energy case), or the LHC (for the high energy case); a much bigger storage 
ring to stack the ions; a longer source to detector distance, to match the first neutrino 
oscillation maximum. New sites for the far detector are to be identified like the Gran 
Sasso Underground Laboratory at about 700 km (for a boost of 350), or the Canfranc 
Underground Laboratory (for 7=150). These scenarios benefit of increased statistics, of 
the energy dependence of the events and eventually of matter effects, with an enlarged 
discovery potential: an increased sensitivity; information on the still unknown neutrino 
mass hierarchy; a reduction of the degeneracies with a better identification of the true 
location of the unknown angle and of the phase. From the physics point of view these 
scenarios with 7 > 150, and particularly the medium energy ones, are very promising. 
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However, they keep, at least for the moment, at a more speculative level than the original 
one, since they are not based on existing accelerator complex, nor on a robust evaluation 
of the ion intensities. Clearly these options deserve further studies. 

The analysis of a CP violating experiment will depend on the knowledge of all 
other neutrino oscillation parameters at the time it takes place. At the moment the 
unknown sign of Am| 3 and octant of # 2 3 introduce supplementary fake regions in the 
multidimensional oscillation parameter space besides those that correspond to the true 
solutions of 6>i3 and 5. A significant reduction of the degeneracies can be obtained in the 
medium and high energy options. On the other hand, the original beta-beam scenario 
has no sensitivity to the octant and to the sign issues, and is affected by degeneracies, 
even combined with a super-beam. The synergy of the accelerator measurements with 
atmospheric data in the same detector, as recently proposed, has been shown to be a very 
promising strategy to get rid of the clone solutions in this case as well. In particular, the 
hierarchy can be identified at 2a CL. from a combined analysis if sin 2 # 13 > 0.03 — 0.05 
for the beta-beam and the super-beam, while the octant can be determined through the 
super-beam, thanks to the atmospheric data. 

The establishment of a beta-beam facility is a natural site to have low energy 
beta-beams with neutrinos in the 100 MeV region. Such options open completely new 
possibilities besides CP, T, and CPT symmetries studies. The physics potential is being 
investigated. It covers different topics from nuclear structure and neutrino-nucleus to 
fundamental interaction studies as well as core-collapse Supernova physics. Since the 
neutrino flux emittance is low, a devoted storage ring appears as more appropriate for 
these applications. 

In conclusion, a beta-beam facility has a rich and broad physics potential. The 
feasibility study of the original scenario, as well as of the small storage for low energy 
beta-beams, is now ongoing within the EURISOL Design Study. Both its outcome and 
future physics studies will furnish the necessary elements to assess the final discovery 
reach. The ongoing work has already shown that a beta-beam facility can clearly offer a 
very competitive strategy for future neutrino accelerator experiments, and much more. 
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